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^~^ Abstract Planetary migration is the process by which a forming planet undergoes a 

'^ drift of its semi-major axis caused by the tidal interaction with its parent protoplan- 

etary disc. One of the key quantities to assess the migration of embedded planets is 
I— I the tidal torque between the disc and planet, which has two components: the Lind- 
Ph blad torque and the corotation torque. We review the latest results on both torque 
M^ components for planets on circular orbits, with a special emphasis on the various 
(— I processes that give rise to additional, large components of the corotation torque, and 
f^ I those contributing to the saturation of this torque. These additional components of 
JL the corotation torque could help address the shortcomings that have recently been 
%^ exposed by models of planet population syntheses. We also review recent results 
j/3 concerning the migration of giant planets that carve gaps in the disc (type II migra- 
te tion) and the migration of sub-giant planets that open partial gaps in massive discs 

(type III migration). 

> 

0^ 1 Introduction 

m 



The extraordinary diversity of extrasolar planetary systems has challenged our un- 
derstanding of how planets form and how their orbits evolve as they form. Among 
^^ the many processes contemplated thus far to account for the observed properties 

T— I of extrasolar planets, the gravitational interaction between planets and their parent 

►jjl protoplanetary disc plays a prominent role. Considered for a long time as the key 

• ^ ingredient in shaping planetary systems, planet-disc interactions, which drive the 
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well-known planetary migration (a drift of a planet's semi-major axis during the 
lifetime of the gaseous disc), have recently been considered by many as being over- 
emphasized. On the one hand, observational data show evidence for vigorous mi- 
gration in many planetary systems, as stressed by the existence of hot Jupiters, Nep- 
tunes and Super-Earths (the recently discovered Kepler-20 planetary system, with 
coplanar rocky and icy planets alternating at periods less than 80 days |30|, provides 
a good example), or by the existence of many mean-motion resonances. Yet, there 
is also compelling evidence that other processes are capable of altering the orbits as 
dramatically as planet-disc interactions (the existence of highly-eccentric or retro- 
grade planets is an example). Also, although many systems seem to have undergone 
orbital migration, many others display planets at distances from their star that are 
of same order of magnitude as the distances of the planets in our Solar System to 
the Sun. One may be tempted to conclude from this that theories of planet-disc in- 
teractions are overrated, to the point that they could be merely ignored in scenarios 
of formation of planetary systems. Yet, as will be detailed in the following sections, 
the following facts are difficult to circumvent: 

• Each component of the torque exerted by the disc on a planet is so large that it 
can half or double the planet's semi-major axis in a time that is usually two or 
more orders of magnitude shorter than the lifetime of protoplanetary discs. 

• These torque components do not cancel out. The residual torque amounts to a 
fair fraction of each torque component, so that one should in general expect that 
planet-disc interactions have a strong effect on planets orbits over the disc life- 
time. 

The central difficulty in planetary migration theories lies precisely in predicting 
the residual torque value. In addition to being a fair difference between several large 
amplitude torques, it is very sensitive to the disc properties near the planet's orbit 
(e.g., density, temperature profiles). This by no means implies that the total torque 
is negligible, but it helps understand why migration theories are slowly maturing. 

One of the main purposes of this review is to provide the reader with an up- 
to-date presentation of the state of planet-disc interactions, with emphasis on the 
torque formulae that govern the migration of low-mass planets. The reasons for this 
special emphasis are three-fold: 

• Low-mass planets are the most critical in planetary population synthesis, as they 
potentially undergo the fastest migration. 

• The sensitivity of detection methods has increased to a point where we can find 
a plethora of Neptune-sized planets or below, which exclusively underwent the 
migration processes typical of low -mass planets during their formation. 

• The subject has been the focus of significant efforts in the recent past. 

The torque acting on a low-mass planet in circular orbit can be decomposed into 
two components: (i) the differential Lindblad torque, arising from material passing 
by the planet at supersonic velocities, which is deflected by the latter and therefore 
exchanges angular momentum and energy with the planet, and (ii) the corotation 
torque, arising from material slowly drifting with respect to the planet, in the vicin- 
ity of its orbit. The differential Lindblad torque has been extensively studied from 
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the early times of planetary migration theories, and is known in much greater de- 
tail than the corotation torque. In fact the corotation region, which has been under 
intense scrutiny over the last five years, has proved to have a much more complex 
dynamics than previously thought. In particular, the value of the corotation torque 
depends sensitively on the radiative properties of the gas disc, and may exhibit large 
values when the gas is radiatively inefficient (as is generally expected in regions 
of planet formation). In addition to this complexity, new problems emerge as the 
computational resources render tractable the task of simulating a planet embedded 
in realistic discs, namely three-dimensional discs invaded by turbulence. 

This review is organised as follows. After a brief description of the physical 
model and notations in § |2] we present in § [5] the migration of low-mass planets 
(type I migration). We detail some recent results on the differential Lindblad torque 
in § |3.1[ and we put special emphasis on the recent developments on the corota- 
tion torque in § 3.2 The migration of gap-opening planets is then examined in §l4j 



with type III migration in § 4.2 followed by type 11 migration in § 4.3 Finally, in 
section §l5l we discuss some recent themes related to planet-disc interactions, such 
as the discovery of massive planets at large orbital separations, and recent models 
of planetary population syntheses. Most sections end with a brief summary of their 
content. 



2 Physical model and notations 

In most of the following we shall consider two-dimensional discs, considering ver- 
tically averaged or vertically integrated quantities where appropriate. At the present 
time, most of the recent investigation on the migration of low-mass planets has been 
undertaken in two dimensions (with a list of exceptions that includes, but is not re- 
stricted to ||T02l|22l|23l[l2l|TIT]|86l[87l|43l|3l), and much insight can be gained into 
the mechanisms of the different components of the torque exerted by the disc on the 
planet through a two-dimensional analysis. It should be remembered, however, that 
two-dimensional results are plagued by the unavoidable use of a softening length 
for the planet's gravitational potential, which fits a two-fold purpose: 

• It mimics the effects of the finite thickness of a true disc, by lowering the magni- 
tude of the planet's potential well. 

• In numerical simulations, it avoids the potential divergence at the scale of the 
mesh zone. 

For this reason, the reader should bear in mind that the ultimate torque expres- 
sions should be sought by means of three-dimensional calculations, and that two- 
dimensional calculations are only used in a first step to elucidate the mechanisms 
that contribute to the torque. At the time of writing this manuscript, most of the 
physics of the torque in two dimensions is fairly well understood, which is why we 
put special emphasis on the two-dimensional analysis. 
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We consider a planet of mass Mp orbiting a star of mass M* with orbital frequency 
Qp. We denote by q the planet-to-primary mass ratio. The planet is assumed to be 
on a prograde circular orbit of semi-major axis a, coplanar with the disc, so that we 
do not consider in this work eccentric, inclined, or retrograde planets. The proto- 
planetary disc in which the planet is embedded is modelled as a two-dimensional 
viscous disc in radial equilibrium, with the centrifugal acceleration and the radial 
acceleration related to the pressure gradient balancing the gravitational acceleration 
due to the central star. We use P to denote the vertically integrated pressure, and s 
to denote (a measure of) the gas entropy, which we express as 

where E is the surface density of the gas and 7 the ratio of specific heats Cp/Cy. 
We denote with T the vertically averaged temperature. We assume in most of what 
follows that the surface density and temperature profiles are power laws of radius, 
with indices a and j3, respectively: 

E oc r" (2) 

and 

Tocr-I^. (3) 

The disc pressure scale length is // = c^/Q, with Q the gas orbital frequency and 
Cs the sound speed. We define the disc aspect ratio hy h^ H/r. When T is a power 
law of radius, so is h, with an index / dubbed the flaring index: 

/10c/, (4) 

which satisfies j3 = 1 — 2/. In almost all studies of planet-disc interactions, the disc 
is modelled with a stationary kinematic viscosity v, aimed at modelling the disc's 
turbulent properties. We will consider that v can be written as v = UyH^Q 11001 , 
with ttv denoting the alpha viscous parameter associated with the turbulent stresses 
in the disc. 

We will also make use of Oort's constants. The first Oort's constant scales with 
the shear of the flow, 

A = -r , (5) 

while the second Oort's constant scales with the vertical component of the vorticity 
of the flow Oj.: 

2r dr 2 ^ ' 

Whenever used, A and B are implicitly meant to be evaluated at the planet's orbital 
radius. 

The governing equations of the flow are the equation of continuity, the Navier- 
Stokes equations and the energy equation (except when dealing with isothermal 
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discs, as specified below), together with the closure relationship provided by the 
equation of state, which is that of an ideal gas. We do not reproduce these governing 
equations here, but refer the interested reader to, for example, ll68l [T4l [84] [T9l l23l 

EH. 



3 Migration of low-mass embedded planets: Type I migration 

Up until recently, type I migration referred to the regime of migration of low-mass 
planets that could be tackled through a linear analysis 146111021 . Recently, however, 
ll89l have shown that one of the torque components, namely the corotation torque, 
can become non-linear at all planetary masses, provided the disc viscosity is suffi- 
ciently small. We shall nevertheless still qualify type I as the migration of low-mass 
planets, up to a threshold mass that we shall specify in § 4. 1 We entertain below the 



two components of the tidal torque: the differential Lindblad torque and the coro- 
tation torque. In the following, we address the properties of type I migration by a 
direct inspection of the tidal torque F, where the planet migration rate a is given by 



IBadMp = r. 



(7) 



3.1 Differential Lindblad torque 



Fig. 1 Disc's surface den- 
sity perturbed by a low- 
mass planet. Streamlines in 
the frame corotating with 
the planet are over plotted 
by solid curves. The set of 
streamlines that librate with 
respect to the planet delimits 
the planet's horseshoe region. 




The differential Lindblad torque accounts for the exchange of angular momentum 
between the planet and the trailing density waves (spiral wakes) that it generates in 
the disc (see illustration in Fig. [Til. The density waves propagating inside the planet's 
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orbit carry away negative angular momentum, and thus exert a positive torque on the 
planet, named the inner Lindblad torque. Similarly, the spiral density waves prop- 
agating outside the planet's orbit carry away positive angular momentum, which 
corresponds to a negative torque on the planet (the outer Lindblad torque). The an- 
gular momentum of a planet on a circular orbit scales with the square root of its 
semi-major axis. The inner Lindblad torque thus tends to make the planet move out- 
wards, while the outer Lindblad torque tends to make it move inwards. The residual 
torque, called differential Lindblad torque, results from a balance between the inner 
and outer torques. In the absence of viscosity, the angular momentum taken away 
by the wakes is conserved until wave breaking occurs, resulting in the formation 
of a shock and the deposition of the wave's energy and angular momentum to the 
disc ll33l . The impact of non-linearities induced by the planet's wakes, which in par- 
ticular lead to the formation of a gap about the planet's orbit, will be described in 
Sectionl4Tn 

The one-sided and differential Lindblad torques can be evaluated in different 
manners: 

1 . In a fully analytic manner upon linearization of the flow equations ||3T1 . In this 
framework, waves propagate away from Lindblad resonances with the planet 
11061, and they constructively interfere into a one-armed spiral pattern Ii83l . 
which begins where the Keplerian flow is supersonic with respect to the planet 
ll33l . One-sided Lindblad torques are then evaluated as the sum of the torques 
arising at each Lindblad resonance. The locations of Lindblad resonances are 
shifted with respect to their nominal location (given by the condition of mean- 
motion resonance with a test particle) by pressure effects. In particular, reso- 
nances with high azimuthal wavenumber have accumulation points at ±2///3 
from the planetary orbit, instead of accumulation at the orbit. This provides a 
torque cutoff Il32l l2l. which can only be evaluated approximately [1 J. This ren- 
ders fully analytic methods of Lindblad torque calculations only approximate. 

2. The torque can be evaluated by solving numerically the linearised equations of 
the flow. This approach was initially undertaken by B6ll . and recently revisited 
by 188] and |84|. 

3. An intermediate approach may be used, in which one solves numerically linear 
equations obtained by an expansion of the flow equations in H /r, where H is the 
disc thickness ||102J . 

Fig.l2]illustrates a number of properties of the Lindblad torque that provide some 
insight into its scaling with the disc and planet parameters. In this figure, the torque 
value has been obtained through the use of analytical formulae similar to that of 
II 1081 (Eqs. 3 to 7 therein), except for the introduction of a softening length for the 
planet's potential, and for the minor correction consisting of the introduction of a 
factor i2 / K" in the forcing terms (see |74|, after their Eq. 13; K" denotes the horizontal 
epicyclic frequency). This figure shows that the torque undergoes a sharp cutoff past 
a peak value, which is found to be of order mmax ^ (2/3) (r///). Also, the dashed 
line shows that, up to the cutoff, the one-sided Lindblad torque approximately scales 
with m^ - as expected from the WKB analysis of |,32J - from which we infer the 
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Fig. 2 Torques at individual outer Lindblad resonances (diamonds) and inner Lindblad resonances 
(triangles), in absolute value. Results are obtained for a disc with aspect ratio H /r = 0.02, and with 
a softening length e = 0.5// for the planet potential. Torques are normalised to the torque value 
given by Eq. (Isj. Stars show the average value of the inner and outer Lindblad torques at a given 
azimuthal wavenumber m, and the dashed line illustrates its m- dependence at small m. 



one-sided Lindblad torque, summed over m, to scale approximately as m^^,^, i. e. as 
{r/H)^ . Besides, the torque naturally scales with the disc's surface density and with 
the square of the planet mass, and dimensional arguments further imply that it ought 
to scale as: 

(8) 



which is indeed the scaling of the one-sided Lindblad torque II108II . 

Fig. |2] also shows that there exists an asymmetry between the outer and inner 
torques, the former being systematically larger than the latter The reasons for this 

. The relative asymmetry is found to scale 



asymmetry are examined in depth in II108I 
with the disc thickness: in thinner discs, the relative asymmetry of the inner and 
outer torque is smaller (which can be understood as due to the accumulation points 
of the inner and outer resonances lying closer to the orbit). As a consequence, the 
differential Lindblad torque scales with: 



The asymmetry between the inner and outer torques also depends upon: 



(9) 



The temperature gradient, since it affects the location of Lindblad resonances, 
and therefore the magnitude of the forcing potential at a given resonance 110611 . 
For instance, a steeper (decreasing) temperature profile decreases the disc's an- 
gular frequency (by increasing the magnitude of the radial pressure gradient), 
which shifts all Lindblad resonances inwards. Outer resonances get closer to the 
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planet orbit, which strengthens the outer torque, whereas inner Lindblad reso- 
nances are shifted away from the planet orbit, which decreases the inner torque. 
The net effect of a steeper temperature gradient is therefore to make the differ- 
ential Lindblad torque a more negative quantity. In the same vein, a shallower 
temperature gradient would shift all Lindblad resonances outwards, making the 
differential Lindblad torque a more positive quantity. The differential Lindblad 
torque may become positive for positive temperature gradients 1 108 1. 

The surface density gradient, which affects the location of Lindblad resonances 
in the exact same way as the temperature gradient, but now also because the 
torque at a given Lindblad resonance directly scales with the underlying surface 



density M106I . A steeper (decreasing) density profile naturally increases the mag- 
nitude of the inner torque compared to that of the outer torque, but this effect is 
mostly compensated for by an inward shift of all Lindblad resonances (just like 
when steepening the temperature profile, as described above). It implies that the 
differential Lindblad torque is quite insensitive to the density gradient near the 
planet location. 

• The disc's self-gravity, which also impacts the location of Lindblad resonances, 
essentially by changing the wave's dispersion relation ||96l[T0ll . 

An accurate determination of the asymmetry between the inner and outer torques 
yields a dimensionless factor to be put in front of i^gf to give the expression for the 
differential Lindblad torque. This issue has triggered a lot of theoretical efforts in 
the last three decades, and it is not completely solved yet. As of the writing of this 
review, the two main results are: 

• An expression obtained by solving numerically the linearised equations of the 
flow in two dimensions with a softened planet potential, in discs with arbitrary 
gradients of surface density and temperature 1881 l84l . It takes the form: 

Fl „. .fOA^'-'' 



I^, (2.5 + 1.7^-O.la)^— ) , (10) 



with a and j3 defined in Eqs. (|2| and (l3|. The expression in Eq. (lOi is most 
accurate for softening lengths e ^ 0.4//. 
• An expression obtained by a semi-analytic method for globally isothermal, three- 
dimensional discs with arbitrary gradients of surface density II 10211 : 

-^=-(2.34-0.10a). (11) 



We note that Eqs. ( 10 1 and ( |TT) exhibit a similar behaviour (for the case ji — 0, 
exclusively contemplated by fl021), that is to say a constant term of similar mag- 
nitude, and a weak dependence on the surface density gradient. Yet, the latter de- 
pends strongly on the softening length, as can be noticed by comparing to the two- 
dimensional, unsmoothed, globally isothermal expression also provided by 11021 : 
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— =-(3.20+ 1.47a). (12) 

This raises the question of whether the dependence on the temperature gradient, 
in a three-dimensional disc, would be as steep as that of Eq. (flOJl. Thus far this 
is an unanswered question, even if recent numerical simulations seem to indicate 
that the dependence of the differential Lindblad torque on the temperature gradient 
in a three-dimensional disc is comparable to that of a two-dimensional disc with 
a smoothing length e ~ 0.4// (Casoli & Masset, in prep). We also comment that 
self-gravity slightly enhances the amplitude of the differential Lindblad torque by 
a factor approximately equal to (1 +2p'), with Qp the Toomre Q-parameter at the 
planet's orbital radius ifTOl . 

Unless the disc has a temperature profile that strongly increases outward, the dif- 
ferential Lindblad torque is a negative quantity which, by itself, would drive type 
I migration on timescales shorter than a few x 10^ yrs for an Earth-mass object in 
a disc with a mass comparable to that of the Minimum Mass Solar Nebula |108j. 



Note however that Eqs. (10 1 and (111 have been derived under the assumption that 
the profiles of surface density and temperature are power laws of the radius. Lo- 
cal variations in the disc's temperature and/or density profiles, due for example to 
opacity transitions fT^ or to dust heating |34| may change the sign and magnitude 
of the differential Lindblad torque. An approximate generalisation of the Lindblad 
torque expression, valid in non-power law discs, has been derived by ll67l for two- 
dimensional discs with a softening length £ = 0.6//. This expression reads: 

"^ -(2.00-0.16a + l.llj3-0.80[i32+-i32"]), (13) 



■T^ef 



where 



h-nrz^^ (14) 



dilogr 

and where a quantity with a ± subscript is to be evaluated in r = azLH/5. Obvi- 
ously, if the temperature profile is a power-law of radius, one has P2 = Pt^ — 0' ^"d 



Eq. ( 13 I reduces to a standard linear combination of a and /3. 

The differential Lindblad torque has also been investigated in strongly magne- 
tised, non-turbulent discs. The case of a two-dimensional disc with a toroidal mag- 
netic field has been studied by 11031 through a linear analysis. She found that the dif- 
ferential Lindblad torque is reduced with respect to non-magnetised discs (as waves 
propagate outside the Lindblad resonances at the magneto-sonic speed (c^ + v^) '' ^, 
with va the Alfven speed, rather than the sound speed c,). Additional angular mo- 
mentum is taken away from the planet by the propagation of slow MHD waves in a 
narrow annulus near magnetic resonances. These results were essentially confirmed 
by ||28l with non-linear 2D MHD simulations in the regime of strong toroidal field 
(the plasma j3 -parameter, /3 — c\jv\, was taken equal to 2 in their study). More re- 
cently, 2D and 3D disc models with a poloidal magnetic field were investigated by 
II79I with a linear analysis in the shearing sheet approximation. While the differen- 
tial Lindblad torque is reduced similarly as with a toroidal magnetic field, extraction 
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of angular momentum by slow MHD and Alfven waves is found to occur in three 
dimensions only. 

We sum up the results presented in this section. 

• The differential Lindblad torque corresponds to the net rate of angular momen- 
tum carried away by density waves (wakes) the planet generates in the disc at 
Lindblad resonances. 

• The sign and magnitude of the differential Lindblad torque arise from a slight 
asymmetry in the perturbed density distribution associated to each wake. 

• The differential Lindblad torque is a stationary quantity, largely independent of 
the disc's turbulent viscosity. Alone, it would drive the migration of Earth mass 
embedded planets in as short a time as a few x 10^ yrs in typical protoplanetary 
discs. 



3.2 Corotation torque 

The other component of the tidal torque, the corotation torque, has long been ne- 
glected in studies of planetary migration. Firstly, it was shown to have a lower 
absolute value than the differential Lindblad torque for typical (decreasing) radial 
profiles of the disc's surface density |102|. Secondly, this torque component, for 



reasons that will be presented at length in section 3.3 should tend to zero after a fi- 
nite time in an inviscid disc (the corotation torque is said to "saturate"). However, as 
indicated by II107I . some amount of turbulence should prevent the corotation torque 
from saturating, and in the last decade the asymptotic value of the torque at large 
time in the presence of dissipative processes has been tackled either analytically 
||54l [5]|, by means of numerical simulations |65|, or both 1691 [85l . Besides, it was 
discovered by |86| that in radiative discs, the corotation torque could, under certain 
circumstances, be so large and positive that it could largely counteract the differ- 
ential Lindblad torque, thereby leading to outward planetary migration. This was 
subsequently interpreted as a new component of the corotation torque arising from 
the disc's entropy gradient ||9l[88ll68ll8l. 

The corotation torque on low-mass planets is usually linked to the so-called 
horseshoe drag, which corresponds to the exchange of angular momentum between 
the planet and its horseshoe region. The planet's horseshoe region encompasses the 
disc region where fluid elements are on horseshoe streamlines with respect to the 
planet orbit (see Fig. [T}. We therefore start by explaining the concept of horseshoe 
dynamics and horseshoe drag. 



3.2.1 Horseshoe dynamics 

In the restricted three-body problem (RTBP), it is useful to write the Hamiltonian 
of the test particle in the frame that corotates with the secondary. In this frame, the 
potential does not depend explicitly on time and the Hamiltonian H is therefore 
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conseiTed. It reads: 



H 



QpL, 



(15) 



where E is the total energy of the test particle as seen in an inertial frame centred 



on the primary, and L its angular momentum. The Hamiltonian of Eq. (15 i is usu 



ally called the Jacobi constant. A sketch of the lines of constant Jacobi value when 
the kinetic energy is zero (E therefore exclusively amounts to the potential energy), 
named zero velocity curves (ZVC), reveals the existence of a horseshoe-like region 
encompassing the secondary's orbit ||78| . Furthermore, in the limit of a small sec- 
ondary's mass, the trajectory of the guiding centre of the test particle in this area is 
shown to have a radial displacement from the orbit that is twice that of the associated 
ZVC ll78l . so that the guiding centre of a test particle can also exhibit a horseshoe- 
like motion in the vicinity of the orbit. In a similar fashion, a gas parcel in a disc 
with pressure can exhibit a horseshoe-like motion, very similar to that of the RTBP. 
There are, however, important differences between these two cases. Firstly, in the 
RTBP, the test particle can exhibit epicyclic motion on top of its global horseshoe- 
like trajectory. In the horseshoe region of a low mass planet, where no shocks are 
present, the gas parcels cannot cross each other's orbits, and therefore essentially 
follow nearly circular streamlines far from the planet. Secondly, the width of the 
horseshoe region is quite different in the two cases. It is much more narrow, for the 
same planetary mass, in the gaseous case than in the RTBP. Also, it has a different 
scaling with the planetary mass in the two cases: in the RTBP, the width of the horse- 
shoe region scales with the cube root of the planetary mass, whereas in the gaseous 
case, provided the planet mass is not too large (we will specify how large below), 
it scales with the square root of the planet mass. To understand the reasons for this 
difference, we depict in the left-hand panel of Fig.[3]the streamlines in the vicinity 
of the coorbital region of a low-mass planet. The planet is located at radius r = 1, 
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Fig. 3 Left: streamlines in the vicinity of an Earth-mass planet embedded in a disc with H /r = 
0.05. The softening length of the potential is e = 0.6//. The thick lines show the separatrices of the 
horseshoe region (the frontiers between the dashed streamlines that exhibit horseshoe-like motion 
and the solid ones corresponding to circulating material with respect to the planet orbit). A few 
arrows show the flow's direction with respect to the planet. The horizontal dashed line shows the 
corotation radius, where the angular velocities of the disc and the planet are equal. Right: effective 
potential at the corotation radius, as a function of azimuth. 
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and azimuth (p — Q. The separatrices of the planet's horseshoe region are depicted by 
thick curves. We see that, quite to the contrary of the RTBP, there is no equivalent to 
the Roche lobe region around the planet (no circulating fluid material bound to it). 
Another difference is that the fixed point (or stagnation point) at which the separa- 
trices intersect lies on the orbit (whereas in the RTBP, they intersect at the Lagrange 
points, away from the orbit, on a line joining the central star and the planet |78|). 
The azimuth of the stagnation point corresponds to the azimuth where, at corotation, 
the effective potential (the sum of the gravitational potential and fluid enthalpy) is 
minimum. This is illustrated in the right panel of Fig. [3] The sign and value of the 
stagnation point's azimuth is closely related to the asymmetry of the inner and outer 
wakes generated by the planet i90| . Note that there may be several stagnation points 
near the planet's corotation radius, depending on the softening length of the planet's 
potential |14|. 

Assuming that the fluid motion is in a steady state in the vicinity of the planet, 
we may use a Bernoulli invariant in the corotating frame ll72ll70l[T4l . This invariant 
can be cast as; 



v2 . 1 



.2r.2 



Bj = - + ri+0--r'n'p, (16) 

where v is the fluid velocity in the corotating frame, rj the fluid's specific enthalpy, 
the sum of the star's and the planet's gravitational potentials, and Qp the planet's 
angular frequency. Equating the value of the Bernoulli invariant at the stagnation 
point (where by definition v = 0), and far from the planet on a separatrix, one finds 
the following expression for the half- width jc^ of the planet's horseshoe region, away 
from the planet: 

x,,oc\cp^, + n'\'J\ (17) 

where 0p is the planetary potential and Tj the perturbation of the gas specific en- 



thalpy introduced by the planet. The s subscript on the right hand side of Eq. ( 17 1 
means that these variables are to be evaluated at the stagnation point. When the 
planet mass is sufficiently small, the streamlines are found to be in good agreement 
with those inferred from the linear expansion of the perturbed velocity field. Note 
that this does not imply that the corotation torque is in general a linear process. 
As was shown indeed by [89 J , the torque exerted by the coorbital material on the 
planet eventually becomes non-linear, no matter how small the planet mass is, pro- 
vided dissipative effects are sufficiently small. However, the fact that important non- 
linear processes take place in this region hardly affects the streamlines themselves. 
In this low-mass regime, the stagnation point has therefore a location independent 
of the planet mass, necessarily at the corotation radius. In that case, the effective 



perturbed potential <Pp + ri scales with the planet mass (Mp), and Eq. ( 17 1 implies 



1/2 

that Xs °<: Mp . This is no longer true when the location of the stagnation point de- 
pends on the planet mass. In particular, in the high-mass regime, the planet gravity 
dominates over the perturbed enthalpy, and the situation resembles that of the RTBP. 
The Bernoulli invariant at the inner and outer stagnation points (which are Li-like 
and L2-like, respectively) is dominated by the planetary potential term and so scales 
as Mp/Rn, where Rh = a{Mp/M^,y'^ is the planet's Hill radius. In the high-mass 
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region, the width of the horseshoe region therefore scales with the cubic root of the 
planet mass, as in the RTBP [70|. 

Lastly, another important difference between the RTBP and the case of a low- 
mass embedded planet is that of the U-turn timescale. Firstly, it should be noted 
that the expression U-turn timescale is ambiguous, for it depends on the horseshoe 
streamline under consideration. The closer to corotation, the longer it takes to per- 
form a horseshoe U-turn, and the U-turn time reaches its minimum value close to 
the separatrix. It is usually this minimum value that is meant by the ambiguous ex- 
pression U-turn timescale. Note that the corotation torque nearly reaches a constant 
value (notwithstanding saturation considerations, that we shall contemplate later) 
after this timescale, as it corresponds to the fluid elements that most contribute to 
the torque, because their angular momentum jump is the largest, and because their 
mass flow-rate is large. 

In the RTBP, the U-turn timescale is of the order of the dynamical timescale Tdyn, 
i.e. a planet orbital period. If one regards, in a crude approximation, the case of a 
low-mass embedded planet as an expurgated version of the RTBP, where most of 
the initial horseshoe streamlines are made circulating, and where only those lying 
close to corotation keep their horseshoe character, one expects the U-turn timescale 
in this case to be significantly longer than the dynamical timescale. This is indeed 
the case: the U-turn timescale is approximately hpT\ib |9j, with hp the disc's aspect 
ratio at the planet's orbital radius, and where T^b is the libration timescale, i.e. the 
time it takes to complete a closed horseshoe trajectory. The libration timescale reads 

Tlib==:j7^ • (18) 

3iipXs 

An alternate, equivalent expression for the U-turn timescale is Tu_tum ^ '^dynH/xs, 
which is corroborated by numerical simulations in which one monitors the advec- 
tion of a passive scalar. This expression shows that the U-turn timescale can indeed 
be longer than the dynamical timescale by a significant factor, as Xs can be much 
smaller than H for deeply embedded, low-mass objects. 
We sum up the results presented in this section: 

• The coorbital region of a low-mass embedded planet in a gaseous disc exhibits a 
horseshoe-like region. 

• This region is much more narrow than in the restricted three-body problem, and 
its radial width scales with the square root of the planetary mass. 

• The stagnation points are located at the corotation radius. There is no equivalent 
to the Roche lobe region for low-mass objects. 

• The horseshoe U-turn timescale is significantly longer than the dynamical timescale. 



3.2.2 Horseshoe drag: an overview 

Far from the horseshoe U-turns in the vicinity of the planet, a fluid element or test 
particle essentially follows a nearly circular orbit, and therefore has a nearly constant 
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angular momentum. When it reaches a U-tum, a fluid element is either sent inward 
or outward, thereby crossing the planet orbit. It does so by exchanging angular mo- 
mentum and energy with the planet. The torque resulting from the interaction of the 
planet with all the fluid elements in the course of performing their horseshoe U-turn 
is called horseshoe drag II107I . 
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Fig. 4 Time evolution of the total torque (sum of the differential Lindblad torque and corotation 
torque) on a Mp = lO^^M^, planet mass embedded in an inviscid isothermal disc. The black curve 
(bottom X-axis) shows the torque evolution over 2000 planet orbits, while the grey curve (top x- 
axis) focuses on the evolution over the first 100 planet orbits. The dashed line shows the value of 
the differential Lindblad torque, and the dotted curve highlights the corotation torque predicted 
with a linear analysis. Taken from L67J . 



Upon insertion of the planet in the disc, it takes some time to establish the horse- 
shoe drag, namely a time of the order of the horseshoe U-turn timescale |89|. This 
is illustrated in Fig. [4] which displays the time evolution of the total torque on a 
Mp = lO^^M,^ mass planet embedded in a thin (h — 0.05) isothermal disc with uni- 
form density profile. The disc is inviscid in this example. Once established (after 
^ 30 planet orbits in our example), the horseshoe drag remains approximately con- 
stant over a longer timescale, which corresponds to the time it takes for a fluid 
element to drift from one end of the horseshoe to the other (that is, about half a li- 



bration time, given by Eq. (18 1). The value of the horseshoe drag that exists between 
the horseshoe U-turn time, and half the horseshoe libration time, is called the un- 
saturated horseshoe drag. Beyond this stage, subsequent U-turns may cause further 
time evolution of the horseshoe drag depending on the disc viscosity, which will be 
described in § 3.3 In the particular case depicted here, where the disc is inviscid, the 



horseshoe drag eventually saturates (it cancels out) after a few libration timescales. 
Until ! 



3.3 we focus on the properties of the fully unsaturated horseshoe drag. 
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3.2.3 Horseshoe drag in barotropic discs 

A hint of the torque exerted by the coorbital material on the planet can be obtained 
by the examination of the perturbed surface density. Nonetheless, this examination 
is rather difficult, because the density perturbation in the planet's coorbital region 
is very small, typically one or two orders of magnitude smaller than the density 
perturbation associated to the wakes. Yet, as can be seen in Fig. [5] an approximate 
subtraction of the wakes density perturbation reveals two regions of opposite signs: a 
region of positive perturbed density ahead of the planet (0 > 0) and a region of neg- 
ative perturbed density behind the planet (0 < 0), which both yield a positive torque 
on the planet. The sign of the perturbed density in the coorbital region depends on 
the background density profile, here it is uniform (a — 0). The largest perturbations 
can be seen to originate near the downstream separatrix in either case (the outer 
separatrix at negative azimuth, and the inner separatrix at positive azimuth), but the 
perturbation is spread radially and extends much beyond the horseshoe region. This 
is to be expected on general grounds: in a barotropic disc, where the gas pressure de- 
pends only on its mass density, any disturbance near corotation excites evanescent 
pressure waves, which extend typically over the disc pressure length scale (here 
H = 0.05a). 




Fig. 5 Perturbed surface density in the coorbital region of an Earth-mass planet in a disc with 
h = 0.05 and a = (unifonn background surface density). The planet is located at r = 1, (p = 0. 
The solid curves show the separatrices of the planet's horseshoe region. In order to remove the 
planet's wakes and to render this map more legible, we have subtracted the density perturbation 
obtained in a situation where no corotation torque is expected (namely a = 3/2, as will be shown 
below). This cancellation is imperfect, however, as the wakes of the two cases are not strictly 
identical. 
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Even if some insight into the corotation torque can be gained by the examination 
of the perturbed density maps, a much more useful quantity is the vortensity (the 
ratio of the vertical component of the vorticity to the surface density, also known 
as potential vorticity), which is materially conserved away from shocks in inviscid, 
barotropic, two-dimensional discs. This is illustrated in Fig.l6] 11071 has evaluated 




Fig. 6 Advection of vortensity in the horseshoe region of an embedded, low-mass protoplanet. As 
in previous panels, the planet is at r = I, <p = 0, and the separatrices of the planet's horseshoe 
region are depicted by white curves. The unperturbed disc's vortensity decreases with radius. The 
(high) vortensity in the inner side of the horseshoe region (at (p < 0, r < 1) is brought to the outer 
side (at <p < 0, r > 1 ) by the outward U-turns behind the planet. Similarly, the (low) vortensity in 
the outer side of the horseshoe region (at ip > 0, r > 1) is brought to the inner side (at (p > 0, r < 1) 
by the inward U-turns ahead of the planet. 



the torque exerted on the planet by test particles embarked on horseshoe motion, by 



making use of the Jacobi invariant of these particles (see section 3.2.1 1. A similar 
calculation can be performed for fluid motion, provided one uses a Bernoulli invari- 
ant by adding the enthalpy to the Jacobi constant fT^^T?!. In both cases one finds 
that the horseshoe drag has the following expression: 



rHs=8|A|B2fl 
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(19) 



where x denotes the radial distance to the planet orbit. In Eq. ( 19 1, the subscript F 



indicates that the (inverse of) the vortensity T. /(O^ has to be evaluated away from 
the planet, in Front of the latter (0 > 0), and the subscript R indicates that it has 
to be evaluated at the Rear of the planet, away from it (0 < 0). The integral in 
Eq. ( [T9] l is usually called the horseshoe drag integral, and in a barotropic disc it can 
be simplified to yield the following expression | 
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Fhs = ^^^y^l4. (20) 



where the quantity 'f, called the (inverse) vortensity gradient for short, is defined 

^^^_Mm, ,2.) 

fllogr 

and can be recast as 3/2 — a for density profiles that can be approximated as power- 
law functions of radius over the radial width of the planet's horseshoe region. In 



Eq. (20 1, all terms are to be evaluated at the planet's orbital radius. This equa- 
tion shows that in two-dimensional barotropic discs, the horseshoe drag cancels 
out when the surface density profile decreases locally as r^^' ^, while it is positive 
for density profiles shallower than r^^'^. For density profiles strongly increasing 
outward, the horseshoe drag can be sufficiently positive to counteract the (negative) 
differential Lindblad torque, and therefore stall the migration of low-mass planets 
ITTTI . Such density jumps may be encountered near the star's magnetospheric cavity, 
or near the inner edge of a dead zone, across which the disc's effective turbulence 
decreases outward (the dead zone refers to the region near the midplane of proto- 
planetary discs that is sandwiched together by partially ionized surface layers). 



The horseshoe drag expression in Eq. (20 1 exclusively holds in the case of 
barotropic discs. Those, naturally, are an idealised concept, and true discs have a 
more complex physics, which yields a more complex expression for the corotation 
torque. However, in any case, as we shall see, a common component of the corota- 
tion torque is given by Eq. (|20]l, so that baroclinic effects yield additional terms to 
this expression. 

We sum up the results presented in this section. In two-dimensional barotropic 
discs, where the gas pressure only depends on the surface density: 

• The horseshoe drag is powered by the advection of the fluid's vortensity along 
horseshoe streamlines inside the planet's horseshoe region. 

• It is proportional to the inverse vortensity gradient across the horseshoe region 
(that is, the quantity 3/2 + d\ogZ / d\ogr for power law discs). It can therefore 
be negative, zero, or positive depending on the surface density gradient across 
the horseshoe region. For typical discs density and temperature profiles, its mag- 
nitude is smaller than that of the (negative) differential Lindblad torque. 



3.2.4 Horseshoe drag in locally isothermal discs 

A long considered framework, both in analytical and numerical studies is that of 
locally isothermal discs, in which the temperature is a fixed function of radius. No 
energy equation is considered in this case, but the flow is no longer barotropic: the 
pressure becomes a function of the density and position (through the temperature). 
The vortensity is no longer materially conserved. Its Lagrangian derivative features 
a source term arising from misaligned density and pressure gradients, or misaligned 
temperature and density gradients lf58l : 
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(22) 



As the temperature gradient has a radial direction and sensibly the same magnitude 
everywhere in the coorbital region, the strength of the source term depends on the 
density gradient: wherever the azimuthal density gradient is large, the source term 
is large. This occurs at the tip of the horseshoe U-turns where we have a strong 
azimuthal density gradient owing to the density enhancement in the planet's imme- 



diate vicinity. The time derivative in Eq. ( 22 1 can be expressed as a derivative with 



respect to the curvilinear abscissa s along the streamline: 



D 
Ds 



vzx vr 



(23) 



where v is the norm of the fluid velocity in the corotating frame, where we have 
used ds = vdt. As the fluid stagnates in the vicinity of the stagnation point (i.e. 
V can be arbitrarily small, provided one chooses a streamline sufficiently close to 



the stagnation point), the source term of Eq. (23 1 formally diverges in the vicinity 
of the stagnation point. The total amount of vortensity created, integrated over the 
horseshoe streamlines, however, remains finite. Fig.|7]shows a vortensity map in the 
vicinity of a low-mass planet for a disc with a = 3/2 (no background vortensity 
gradient) and j3 = 1 (uniform aspect ratio). 
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Fig. 7 Vortensity field in the coorbital region of a low-mass planet (left) and radial profile of 
vortensity at (|) = +0.5 rad (right), 30 orbital periods after the insertion of the planet. The disc has 
no background vortensity gradient, and a flat aspect ratio H /r — 0.05 . Thin stripes of vortensity 
of opposite signs are clearly visible at the downstream separatrices. One can also see a mild pro- 
duction of vortensity in the wake, which fades away as one recedes from corotation, because of the 
winding of the wake and because fluid elements move faster away from corotation. Contrary to the 
adiabatic case that we shall present in § |3.2.5| the vortensity cut is not singular at the separatrix, as 
can be seen in the right panel. The radial resolution in the run presented is 9.3 x 10^'a. 



Thus far there is no rigorous mathematical proof that the horseshoe drag expres- 



sion of Eq. ( 19 1 still holds in locally isothermal discs, as all demonstrations of that 



expression rely on the existence of a Bernoulli invariant, which does not exist in 
the locally isothermal case. Yet, data from numerical simulations suggest that this 
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expression is still valid in that regime. Assuming its validity from now on, we in- 
fer that the horseshoe drag must exhibit a dependence on the temperature gradient. 
The rational for this being that the outgoing vortensity, accounted for by the horse- 
shoe drag integral, includes the vortensity produced in the vicinity of the planet, 
which depends on the temperature gradient. Two-dimensional numerical simula- 
tions have confirmed the existence of an additional component of the corotation 
torque that depends on the temperature gradient [6 14 84 1. The sign and value of 
this temperature-related corotation torque have a complex, and rather steep depen- 
dence with the softening length of the planet's potential fT4'|. Indeed, the topology 
of the horseshoe region depends heavily on the softening length: at large softening 
lengths, only one X-stagnation point is observed at corotation, in the planet vicinity, 
whereas at low softening lengths, two X-stagnation points are usually observed at 
corotation, on either side of the planet |14, 90 1. The vortensity produced along a 
streamline depends on the path followed by the streamline. This situation is there- 
fore much more complex than in the barotropic case or the adiabatic case that we 
present below, in which the existence of invariants under certain circumstances al- 
lows to get rid of the dependence on the actual path followed by fluid elements 
during their U-turns. Quite ironically, the locally isothermal case, which has served 
as a standard framework for more than two decades, is very difficult to tackle ana- 
lytically. 

The steep dependence of the temperature related corotation torque on the soften- 
ing length appeals for a three-dimensional study of this torque, which is not plagued 
by softening issues. Such study has been undertaken by Casoli & Masset (in prep.), 
who find a linear dependence of the three-dimensional horseshoe drag on the tem- 
perature gradient j5 = —dlogT /dlog r, as steep as the dependence on the vortensity 
gradient Y given by Eq. ([20|. 



3.2.5 Horseshoe drag in adiabatic discs 

In the previous sections, the set of governing equations of the fluid did not include 
an energy equation, and the disc temperature, which was set as a prescribed function 
of radius, did not evolve in time. The first calculations undertaken with an energy 
equation were those of iTTTll and ||86l . The former were devised in the shearing sheet 
framework, so that no net torque could be experienced by the planet, owing to the 
symmetry properties of the shearing sheet. Still, these authors found that radiative 
cooling could significantly affect the perturbed surface density pattern associated 
with the wakes, thus changing the magnitude of the one-sided Lindblad torque. Il86l 
considered a planet in a three-dimensional disc, with an energy equation and ther- 
mal diffusion, and nested grids around the planet to achieve a very high resolution. 
They found that the migration of a low -mass planet could be reversed in sufficiently 
opaque discs, under the action of the corotation torque. The same result was in 
particular obtained in the adiabatic limit, which we are now going to focus on, as 
thermal diffusion yields an additional complexity, not needed at this stage. We will 



take thermal diffusion and other dissipative processes into account in § 3.3 
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It was soon realised that the results of |86| were due to a new component of 
the corotation torque, linked to the entropy gradient (l9][88l. This is illustrated in 
Fig. IS] in which we compare the torque results for 69 different random disc profiles 
(the density slope a being a random variable uniformly distributed over the interval 
[^3/2, +3/2], and the temperature slope j3 being an independent random variable 
uniformly distributed over the interval [—2, +2].) Each calculation has a smoothing 
length of the planet's potential e ~ Q.3H, and an aspect ratio at the planet location 
hp = 0.05. 

For each pair of a and j3, we ran two calculations: a locally isothermal one, and 
an adiabatic one with a ratio of specific heats y— 1 .4. The torque difference, dubbed 
adiabatic torque excess, is then obtained by: 

^FnT-r^d-^. (24) 

The correction of the isothermal torque fiso by a factor y^ ' is necessary as both 
the differential Lindblad torque and the barotropic part of the horseshoe drag (the 
vortensity-related corotation torque) scale with the inverse square of the sound 
speed, which turns out in the adiabatic case to be df^ = cf°y^'^. The right part 
of Fig. IS] shows a clear one-to-one relationship between the adiabatic torque excess 
and the entropy gradient, irrespective of the individual values of a and j3, which 
justifies the name entropy -related corotation torque given to the adiabatic torque 
excess. This torque is shown in two different regimes: the linear regime, soon after 
the insertion of the planet in the disc, and the horseshoe drag regime, reached after a 
longer timescale, as discussed in previous sections. As pointed out by lISSJI . the non- 
linear corotation torque can be much larger (in this example, by a factor of about 5) 
than its linear counterpart, depending on the planet's softening length. 

An interpretation of the entropy-related corotation torque was given soon after 
its discovery by |9l and fSSl. Although this early interpretation was shown not be 
quite correct when the horseshoe drag expression in the adiabatic case was worked 
out subsequently, we describe it here briefly as the mechanism on which it is based 
allows for a clear understanding of the correct origin of the entropy torque. 

In an adiabatic disc, the entropy is materially conserved along the path of fluid 
elements, as long as they do not cross shocks, in a strict analogy with vortensity 
for barotropic flows. If we assume, for instance, a disc that has a positive radial 
entropy gradient prior to the planet insertion, outward horseshoe U-turns (behind the 
planet) bring to the outer side of the horseshoe region the low entropy of the inside. 
Similarly, inward horseshoe U-turns (ahead of the planet) bring to the inner side of 
the horseshoe region the high entropy of the outside. As the disc maintains a pressure 
equilibrium, the relative variations of the pressure across the coorbital region can be 
neglected, so that the surface density features relative variations opposite to that of 
the entropy, by virtue of the first order expansion of Eq. ([T}: 

r^ + -^^-0. (25) 
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Fig. 8 Torque difference between the adiabatic and locally isothermal cases, given by Eq. |24|, as 
a function of the entropy gradient (right panel), at early times (linear stage, triangles) and during 
the horseshoe drag stage (diamonds). The quantity ,y in x-axis is defined in Eq. pOl. Results 
have been obtained with calculations with random values of the surface density slope {—a) and 
temperature slope (— /3), shown in the left plot. The total number of runs is 138. 



As a consequence, two lobes of perturbed surface density appear in the horseshoe 
region ID [88] |44] |43| . that both yield a torque of same sign. These two lobes are 
shown in Fig. |9] 
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Fig. 9 Surface density perturbation in the planet's coorbital region of an adiabatic disc (left) and 
perturbation of entropy (right). The sign of the density lobes is opposite that of the entropy lobes. 
In addition to the lobes, the left plot also shows the wake, saturated in this representation. The 
background entropy profile increases with radius in this example. 



The early interpretations identified the entropy related torque as the torque aris- 
ing from the above density lobes. This explanation was appealing at first, because it 
gives the correct sign for the entropy related corotation torque, and the correct or- 
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der of magnitude: Q performed an approximate, horseshoe-drag like integration to 
evaluate the impact of these lobes on the torque, whereas |[88l performed an approx- 
imate direct summation. Both results were in rough agreement with the magnitude 
of the adiabatic torque excess, at least for the value of the potential softening length 
used in these studies. Yet, this explanation of the entropy torque quickly turned out 
not to be fully satisfactory, for the following reasons: 

• Numerical explorations performed at different smoothing values showed that the 
entropy torque was approximately scaling as e^', down to very low values of the 
smoothing length (e ~' 0.05//). This apparent divergence of the entropy torque 
at low smoothing was incompatible with a scaling of the torque in x^, since the 
half-width x^ of the horseshoe region remains finite in the limit of a vanishing 
softening length 1 90 1 . 

• The saturation of the entropy torque was also problematic. We will examine sat- 



uration processes in detail in § 3.3 but for our purpose it suffices to know that 
the corotation torque always saturates in inviscid discs, that is to say tends to 
zero after a few libration timescales (as we have seen in Fig.|4]l. Therefore, one 
could devise a setup with an inviscid disc and finite thermal or entropy diffusion, 
which would forever maintain the same entropy perturbation within the horse- 
shoe region (finite thermal diffusion is required in order to avoid phase mixing of 



entropy, as we shall see in § 3.3.2 1. As expected, the entropy related corotation 
torque is found to saturate as the disc is inviscid, while the (approximate) same 
density lobe structure is maintained within the horseshoe region [69 1. This im- 
plies the density lobe structure would exert a torque at early times, but not at late 
times, which is contradictory. 



• Finally, as we have seen in § 3.2.3 the density perturbation responsible for 
the corotation torque is not bound to the horseshoe region, but can extend fur- 
ther radially by the excitation of evanescent waves. In the barotropic case, the 
vortensity-related horseshoe drag is actually fully accounted for by an evanescent 
density distribution within the coorbital region. In the adiabatic case, attributing 
the entropy torque {i.e. the whole difference between an adiabatic and an isother- 
mal calculation) to the density lobes was therefore tantamount to assuming that 
the evanescent wave structure in the coorbital region was the same in the adia- 
batic and isothermal cases, which is not obvious. 

The identification of a convenient invariant of the flow for adiabatic discs with 
uniform temperature profile allowed f68l l69J to demonstrate that the horseshoe drag 



expression was exactly the same as that of the barotropic case, given by Eq. ( 19 1. 
In this case, the evaluation of the horseshoe drag amounts again to a budget of the 
vortensity entering or leaving the vicinity of the planet on horseshoe streamlines. An 
important consequence of this is that the torque due to the density lobes must not 
be incorporated manually, separately, into the corotation torque expression, and the 
whole problem of determining the horseshoe drag amounts to an evaluation of the 
vortensity distribution within the horseshoe region. Before we clarify this point, we 
stress that the vortensity distribution within the horseshoe region has the following 
features: 
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Since the flow is baroclinic, vortensity is not materiafly conserved along stream- 
lines. However, contrary to the locally isothermal case, the existence of a flow 
invariant in adiabatic discs with flat temperature profile allows an estimate of the 
vortensity acquired by a streamline during a U-turn, independently of its actual 

pathllMl- 

The vortensity created over the interior of the horseshoe region is very small, and 

has no impact on the torque, because it has same sign on both sides of the planet 



• The main difference arises from a (formally) singular production of vortensity 
(or vorticity) on downstream separatrices, due to the entropy discontinuity at this 
location (which results from the entropy advection within the horseshoe region). 
This (formally) singular production of vortensity is readily apparent in the source 



term of Eq. (22i, and is illustrated in Fig. 10 It can be evaluated analytically 
either using the flow invariant introduced in |68|, or directly using Eq. (J22]i as 
in 11841 . The first approach is self-contained and yields the amount of singular 
vortensity as a function of the flow properties at the stagnation point. The second 
is not restricted to flat temperature profiles, but it requires knowledge of the fluid 
velocity along horseshoe streamlines, which depends on the exact geometry of 
the horseshoe region, much like in locally isothermal discs discussed in § 3.2.4 
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Fig. 10 Vortensity map in the coorbital region of an Earth-mass planet embedded in an adiabatic 
disc with H = 0.05 at the orbital radius of the planet (left), and radial profile of vortensity at (j> = 
0.5 rad. (right), 60 orbital periods after the insertion of the planet in the disc. The grid resolution 
and disc gradients are the same as in Fig.p](a = 3/2, /3 = 0). The vortensity peak has a much more 
compact profile than that of Fig.^ This is to be expected as we have a singular vortensity sheet in 
the adiabatic case, and a continuous one, peaked at the separatrix, in the locally isothermal case. 



We now clarify the contribution of the density lobes to the corotation torque. The 
corotation torque is directly related to the density perturbation within the corotation 
region, which can be written as 
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where, counter to Eq. (25 1, we shall not assume 5P = 0. The pressure perturbation 

5P can be shown to satisfy a second-order partial differential equation ||68l . with 

solution of the form: 

5P 8u 

— ^7K^—, (27) 

P u 

where u = s'" x E/co,, K ^ exp(— j.^cj///) is a Green kernel normalised to unity 
(jK(x)dx — 1, with X the radial distance to the planet orbit), and H is the local 



pressure scale height. Further denoting the inverse vortensity by / = E/co-, Eqs. ( 26 1 



and p7| ) yield 
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The first term on the right-hand side of Eq. ( 28 1 corresponds to the perturbed surface 
density associated to evanescent pressure waves (like in the barotropic case, where it 
reduces ioK* 51 /I), and the second term to the density lobes resulting from entropy 



advection. Since the convolution by the unitary function K in Eq. ( 28 i does not 



change the linear mass of the perturbation (/ 5E{x)dx), the corotation torque is the 



same as if, in the expression for the density perturbation in Eq. (28 1, the convolution 



product were actually discarded 11681 . In the barotropic case for instance, this leads 



to the torque expression given by Eq. ( 19 1. In the adiabatic case, it shows that. 



counter-intuitively, the density lobes exert no net corotation torque. This further 
explains why, akin to the barotropic case, the calculation of the corotation torque 
comes to evaluating the vortensity distribution within the horse shoe region. Since 
the main difference in the vortensity field between the adiabatic and barotropic cases 
is the appearance of a singular sheet of vorticity at the downstream separatrices, and 
given that the magnitude of this sheet scales with the entropy gradient, this singular 
vorticity sheet can be unambiguously identified as the origin of the entropy -related 
torque. 

Upon evaluation of the magnitude of the vorticity sheet at the separatrices, ||68l 
inferred the following expression for the entropy -related corotation torque: 

^FhT = - l^znlq^'h-^ (29) 

where the above expression has been derived in the framework of a flat temperature 



profile, and assuming a ratio of specific heats 7 = 1 .4. In Eq. ( 29 1, all disc quantities 



are to be evaluated at the planet's orbital radius, and the quantity ,5^ is defined by 

=^=-1^' (30) 

7alogr 

and can be recast as [j3 + (7— l)a)]/7for surface density and temperature profiles 
that can be approximated as power-law functions of radius over the planet's horse- 
shoe region. 

Considering discs with arbitrary temperature profiles, f84l also evaluated the pro- 
duction of vortensity at downstream separatrices, which required estimating the ve- 
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locity along streamlines through a fit of numerical simulations. Unlike |68|, they 
manually added the torque contribution from the density lobe structure. In the end, 
the latter remains small compared to the torque contribution from the singular sheet 
of vorticity. This explains why, overall, the derivations of the entropy-related coro- 
tation torque by ll68l and ll84i l are in broad numerical agreement, within 30 %. 
The generalisation to an arbitrary temperature profile of Eq. ([29]l cannot be tack- 



led fully analytically, much as in the locally isothermal case. Yet, Eq. (24i shows 



that the adiabatic corotation torque is the sum of the entropy related term, given by 



Eq. (29 1, and the locally isothermal corotation torque (corrected by a factor 7). The 



latter is itself made up of two terms, as we have seen in sections 3.2.3 and 3.2.4 The 



corotation torque is therefore, in a general situation, the sum of three terms: 

• The vortensity related torque, proportional to the vortensity gradient, and given 
by Eq. ^. 

• The temperature related torque, proportional to the temperature gradient, dis- 
cussed in section [3231 



• The entropy related torque, proportional to the entropy gradient, given by Eq. ( 29 1 



It can be observed that there are only two degrees of freedom for the disc pro- 
files (the density and temperature gradients a and j3 , or the vortensity and entropy 
gradients f and =5^, etc.), so that for a specific setup one may simplify the torque 
expression as a linear combination of the two independent parameters. This simpli- 
fication is not desirable, however, because it blurs the distinct physical origin and 
characteristic of each of the three terms. Besides, one can disentangle these three 
terms by varying parameters such as the smoothing length e or the ratio of specific 
heats 7. Any simplification of the torque expression is thus highly setup dependent. 

Finally, we summarise the main message to take away about the corotation 
torque. 

• In all cases it features a term that scales with the gradient of vortensity across the 



horseshoe region, given by Eq. ( 20 1. It has one or two additional terms, depending 
on whether an energy equation is taken into consideration. The first of those 
scales with the temperature gradient, and if an energy equation is included, there 
is a second one that scales with the entropy gradient. 

In all cases the corotation torque comes from the vortensity distribution in the 
horseshoe region. The additional contributions arise from the vortensity created 
by the temperature gradient and/or the entropy gradient. 



3.3 Saturation properties of the horseshoe drag 



We have described in Section |X2| the physical origin and properties of the corotation 
torque in inviscid discs, with a special emphasis on the fully unsaturated horseshoe 
drag, which is the maximum value the corotation torque may take. This value is 
obtained about one horseshoe U-turn timescale after the planet insertion in the disc, 
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and is maintained over about half a libration timescale. Its sign and magnitude are 
determined by the gradients of vortensity, temperature and entropy across the horse- 
shoe region. 

In the absence of diffusion processes, after about half a libration timescale, the 
vortensity and entropy advected along the downstream separatrices of the horseshoe 
region reach the planet again, undergo another U-turn, and phase mixing starts to oc- 
cur. Vortensity and entropy are progressively stirred up within the horseshoe region, 
and the horseshoe drag oscillates with time with a decreasing amplitude, as shown 
in Fig.ffl The horseshoe drag ultimately cancels out as both vortensity and entropy 
get uniformly distributed after several libration times 1851 l69l . This is known as the 
horseshoe drag saturation. 

Diffusion processes (viscosity, thermal diffusion) may maintain respectively the 
vortensity and entropy gradients across the horseshoe region, and thus sustain the 
horseshoe drag to a non-vanishing value. We review below the saturation properties 



of the horseshoe drag in barotropic discs (§ 3.3. 1 1 and in radiative discs (§ 3.3.2 1. 



3.3.1 Saturation properties of the vortensity-related horseshoe drag in 
barotropic viscous discs 

In barotropic discs, the horseshoe drag saturates as vortensity is strictly advected 
along horseshoe streamlines. Viscosity acting as a diffusion source term in the 
vortensity equation can sustain a non-zero vortensity gradient across the horse- 
shoe region. The vortensity-related horseshoe drag then attains a steady-state value, 
which arises from a net exchange of angular momentum between the horseshoe re- 
gion and the rest of the disc 1641 l65l . This steady-state value depends on how the 
viscous diffusion timescale across the horseshoe region (Tyisc) compares with the 
horseshoe libration timescale (Tub) and the horseshoe U-turn timescale (Tu-tum)- 
Denoting by Vp the kinematic viscosity at the planet location, Tyisc ^ ^1/^p- The 



libration timescale is given by Eq. ( 18 1, and the U-tum timescale is typically a frac- 
tion H /r of the libration timescale |9 1. 

For the corotation torque to remain close to its maximum, fully unsaturated value 
in the long term, the inequality 

Tu-turn < Tyisc < Tlib/2 (31) 

should be verified. When the second inequality is satisfied, the vortensity at the up- 
stream separatrices is kept stationary, which prevents phase mixing of vortensity 
within the horseshoe region |5 , 64 1. When the first inequality is satisfied, vortensity 
is approximately conserved along U-turns, which maximises the effective vorten- 
sity gradient across the horseshoe region f65l l69l l85l . Taking Xj ^ lAa^Jq/hp (as 



measured with a planet softening length w 0.6H), inequality ( 31 1 may be cast as 



0.32^^/2/^^^7/2 ^ ^^^ ^ 0.16^^/2/,^ 9/2^ ^22) 
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where tty.p and hp denote the disc's alpha viscosity and aspect ratio at the planet 
location, respectively. The alpha viscosity for which the corotation torque takes its 
maximum value can be approximated as 0.16^^'^/ip^ |8j. 
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Fig. 11 Top: steady-state torque on a Mp = 6 x lO^^M^, planet mass embedded in a thin (h = 0.05) 
disc for various alpha viscous parameters at the planet location. In this series of runs, a = 1 /2 and 
j8 = 0. Different saturation regimes of the corotation torque are illustrated, depending on how the 
viscous timescale across the planet's horseshoe region (Tyisc) compares with the horseshoe U-turn 
timescale (Tu-tum) and the libration timescale (Tub). The final torque value in an inviscid case, 
which reduces to the differential Lindblad torque, is shown by a dashed line. The fully unsaturated 
total torque (differential Lindblad torque plus fully unsaturated horseshoe drag) in an inviscid run 
is depicted by a dotted line. Bottom: vortensity distribution inside the planet's horseshoe region for 
the three alpha viscosities shown by squares in the top panel (viscosity increases from left to right). 
The separatrices of the horseshoe region are depicted by solid curves, and the planet position by a 
filled circle. 



The saturation properties of the corotation torque are illustrated in Fig. 1 1 for a 
2 Earth-mass planet embedded in a thin (hp — 0.05) viscous disc. The background 
temperature profile is uniform, and the surface density decreases as r^^'^. The top 
panel displays the steady-state torque at different alpha viscosities (a constant kine- 
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matic viscosity V was used in the simulations). The left-hand term in inequality ( [32| l 
is « 1.7 X 10^^, while the right-hand term in sa 1.7 x 10^^, and it is clear from 



Fig. 1 1 that the corotation torque is maximum between these two alpha viscosi- 
ties. When the viscosity is small enough so that Tyisc ^ ^ub, viscosity is ineffi- 
cient at restoring the vortensity gradient across the horseshoe region, and the horse- 
shoe drag takes very small values (it saturates). At very large viscosities, such that 
Tvisc ^ Tu-tum, the corotation torque plateaus at its value in the linear regime ll89ll . 
The vortensity distribution inside the horseshoe region for each saturation regime 
is shown in the bottom panels of Fig. Ill] In the left panel, tty.p ^ 10^^, and the 
steady-state vortensity distribution within the horseshoe region is uniform, result- 
ing in a vanishing corotation torque. In the middle panel, tty.p ^ lO^'* maintains a 
maximum vortensity contrast between the rear and front parts of the planet, with 
the consequence that the corotation torque is close to its fully unsaturated value. 
In the right panel, av,p ^ 10^^ imposes the initial (unperturbed) vortensity profile 
along the horseshoe U-turns, and the horseshoe drag therefore reduces to the linear 
corotation torque. 

All attempts to capture analytically the saturation of the corotation torque have 
been carried out using a simplified streamline model that assumes the drift of the 
coorbital material with the velocity of the unperturbed disc, and which does not re- 
solve spatially nor temporally the U-turns. This model was proposed in |64| and in 
a more formal manner in |69|, where a numerical implementation of it is also de- 
scribed. All analytic works on the corotation torque saturation, whether they provide 
an asymptotic torque value |64 69, 85 1 or they capture the time dependence of the 
torque in an inviscid disc 1 109] make use of this simplified model. Solving for the 
torque asymptotic value in this simplified advection-diffusion model can be tacked 
in a variety of ways. One such way consists of neglecting the azimuthal variation of 
the vortensity so as to reduce the advection-diffusion problem essentially to a one 
dimensional radial problem. This is the approach of |64| and |85|. These two works 
are quite different in their assumptions, and suffer from orthogonal restrictions: 

• 164 1 exclusively contemplates the case of a disc with flat profiles of surface den- 
sity and kinematic viscosity, so that his results must be rescaled by hand to apply 
to a general case. The approach used in this work considers the global angular 
momentum budget of the trapped horseshoe region, and relies upon the eval- 
uation of the viscous friction of the disc on the separatrices. It also takes into 
account the viscous drift of material across the horseshoe region. 

• 185 1 consider a disc with an arbitrary surface density gradient, and directly use 



the horseshoe drag integral of Eq. ( 19 1. Their model assumes no radial drift of 



disc material across the horseshoe region. 

Quite remarkably, these two approaches yield the exact same result, which can be 
cast either in terms of Airy functions |64| or in terms of Bessel functions Il85l . 

One can also solve the advection-diffusion equation satisfied by the fluid's 
vortensity in two dimensions, the solution being exact in the limit of a small vis- 
cosity. In this limit, the problem amounts to an alternation of convolutions (viscous 
diffusion of vortensity between two successive horseshoe U-turns) and reflections 
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(mapping of vortensity - or, vortensity conservation - from one tip of the horseshoe 
region to the other during a U-turn). This is the approach of 1691 . who also discard 
the possible radial drift of disc material across the horseshoe region. The depen- 
dence thus obtained — equation (119) of |69 1 — is broadly the same as that of ll64l 
and ll85l . but reproduces more closely the results from numerical simulations. We 
note that the decay of the torque value found at large viscosity (see Fig. [TT] i, which 
corresponds to a decay towards the linear corotation torque value [89], has not yet 
been described analytically in a self-contained manner. |69| and 1851 use an ad-hoc 
reduction factor, either with one free parameter |69 1 or two free parameters ||85J . the 
value of the free parameters being inferred from numerical simulations. 

In summary, in barotropic discs, vortensity essentially obeys an advection- 
diffusion equation in the coorbital region. When the viscous diffusion timescale 
across the horseshoe region is: 

• Long compared to the libration period, vortensity is progressively stirred up and 
the corotation torque ultimately saturates (tends to zero). 

• Short compared to the libration period, but long compared to the horseshoe U- 
turn time, the corotation torque is close to its fully unsaturated, maximum value. 

• Short compared to the horseshoe U-turn time, the corotation torque tends to its 
value predicted in the linear regime. 



3.3.2 Saturation properties of the horseshoe drag in radiative discs 

Much as in barotropic discs, the estimate of the asymptotic corotation torque value 
in radiative discs amounts to the determination of the vortensity distribution within 
the horseshoe region at later times. There is an additional complexity, however, due 
to the fact that this is no longer an advection-diffusion problem, but an advection- 
diffusion-creation problem, as vortensity is created during the U-turns (see sec- 
tion [ 



3.2.5 and in particular Fig. 10 1. Furthermore, the amount of vortensity created 



depends on the entropy distribution, as was explained in section 3.2.5 This analysis 
was undertaken by |85| in the case of a unitary thermal Prandtl number (the vis- 
cosity V and thermal diffusion % have same value). A corotation torque expression 
was proposed by these authors, as a result of a fit of numerical simulations. Under 
the assumption of a unitary Prandtl number, the parameter space to be explored is 
one-dimensional, and for a (common) value of v and %, the radiative torque is found 
to saturate more easily than the barotropic torque. This is interpreted by the authors 
as due to the fact the entropy-related corotation torque is essentially due to a unique 
streamline, where the advection speed is maximal (that of the separatrices). 

To relax the assumption of a unitary Prandtl number, one may assume that the 
torque dependence upon viscosity or thermal diffusion have the same shape, which 
allows to propose a formula with two independent parameters v and %, which can 
then be validated by checking its accuracy with numerical simulations. This is the 



approach adopted by ||851. Another solution consists in using a streamline model 



such as the one outlined in section 3.3. 1 This is the approach of 1691 . As the vorten 



sity is now determined by an advection-diffusion-creation problem, one needs to 
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amend the barotropic model of section |3.3.1| by adding the creation of vortensity 
during the U-turns, which is determined by the entropy field. Therefore, prior to the 
determination of the vortensity distribution, an analysis of the entropy distribution at 
later times is required. This preliminary determination can be done easily, because 
the entropy obeys an advection-diffusion problem formally equivalent to the vorten- 
sity distribution in the barotropic case, in which one replaces the vortensity with the 
entropy, and the viscosity v with the thermal diffusion %. Once the entropy distribu- 
tion within the horseshoe region is known, the vortensity distribution at late times is 



obtained, which allows, upon the use of the horseshoe drag expression of Eq. ( 19 1, 
for an expression of the corotation torque as a function of viscosity and thermal 
diffusivity, and which can be checked a posteriori against numerical simulations. 

The corotation torque expressions, as a function of viscosity and thermal diffu- 
sivity, are given by Eqs. (161-164) of 169J, or by Eqs. (52-53) of HH. 



3.4 Type I migration in turbulent discs 

We have examined in the previous sections the properties of planet-disc interactions 
assuming viscous discs, described with a stationary kinematic viscosity aimed at 
modelling their turbulent transport properties. Because the corotation torque may 
play a dramatic role in the orbital evolution of low-mass planets, and its magnitude 
is intimately related to diffusion processes taking place within the planet's horseshoe 
region, it is relevant to determine how turbulence may impact type I migration. 

Turbulence in protoplanetary discs can have a variety of origins. These include 
hydrodynamic instabilities, such as Rossby-wave instabilities |58|, the global baro- 
clinic instability Ii421 l59l. the sub-critical baroclinic instability ||49il, planetary gap 
instabilities ll55l l56ll (which we will discuss in § |4.2| i, or the Kelvin-Helmholtz in- 
stability triggered by the vertical shear of the gas as dust settles into the mid plane 
lf41J . Convective instability might also be relevant in the inner parts of massive discs, 
and it would be interesting to examine its impact on type I migration. Perhaps the 
most likely source of turbulence in protoplanetary discs is the magnetohydrody- 
namic (MHD) turbulence resulting from the magnetorotational instability (MRI) 
H. It relies on the coupling of the ionised gas to the weak magnetic field in the disc, 
lonisation may occur in the vicinity of the central object due to the star's irradiation, 
or further out in the disc layers, most probably through the UV background or cos- 
mic rays. It is currently debated which regions of planet formation near the disc mid 
plane are sufficiently ionised ('active') to trigger the MRI, and which ones remain 
neutral (which is usually referred to a 'dead zone'). In the latter case, some trans- 
port of angular momentum would still be present through the propagation of waves 
induced by MHD turbulence in the disc's upper layers |27|. The alpha viscous pa- 
rameter associated to MHD turbulence is typically in the range [5 x 10^^—5 x 10^^] 
in active regions, while being about two orders of magnitude smaller in dead zones. 

The properties of type I migration in weakly magnetised turbulent discs have 
been investigated in a couple of studies. 1, 82 J performed 3D simulations of locally 
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isothermal discs fully invaded by MHD turbulence. They found that the running 
time-averaged torque on a fixed protoplanet experiences rather large-amplitude os- 
cillations over the reduced temporal range over which the simulation could be run, 
and that its final value differs quite substantially from the torque value expected in 
viscous disc models. Similar results were obtained by |[80l . who allowed the planet 
orbit to evolve. A primary reason for the observed difference between the viscous 
torque and the time-averaged turbulent torque is that the 3D MHD simulations were 
not converged in time. This was suggested by ||8l, who considered 2D isothermal 
discs subject to stochastic forcing, using the turbulence model originally developed 
by |48|. They showed indeed that when time-averaged over a sufficiently long time 
period, which may be as long as a thousand orbits, both the differential Lindblad 
torque and the corotation torque behave very similarly as in equivalent viscous disc 
models. These results were essentially confirmed by the 3D MHD simulations by 
Q, who adopted a locally isothermal disc model with a mean toroidal magnetic 
field, in which non-ideal MHD effects and vertical stratification were neglected 
(see illustration in Fig. 12 1. Similar agreement was obtained by |104| with verti- 



cal stratification. Nonetheless, [71 found an additional corotation torque with mod- 
erate magnitude in their 3D MHD simulations, related to the presence of a mean 
toroidal magnetic field. The existence and properties of this additional corotation 
torque have been explored by Guilet, Baruteau & Papaloizou (in prep.) in 2D weakly 
magnetised, non-turbulent disc models, in which the effects of turbulence are mod- 
elled by viscous and magnetic diffusivities. They find that the additional corotation 
torque can take large values, and even exceed the differential Lindblad torque, de- 
pending on the disc's viscous and magnetic diffusivities, and the amplitude of the 
background magnetic field. 




Fig. 12 Perturbation of the disc's density by an embedded planet. Left: adiabatic two-dimensional 
disc. The density lobes within the coorbital region, which arise from the advection of entropy, help 
identify the (tiny) radial width of the planet's horseshoe region. Middle: case of an isothermal 2D 
disc with turbulence induced by stochastic stirring. Right: case of an isothermal 3D disc invaded by 
MHD turbulence due to the MRI (the density in the disc mid plane is displayed). In the middle and 
right panels, the turbulent density perturbation is comparable to the perturbed density associated to 
the planet's wakes. Images taken from I6|&i7j, respectively. 
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The aforementioned results were for embedded planets with a horseshoe radial 
width that is a moderate fraction of the disc's pressure scaleheight, the latter being 
the typical size of turbulent eddies. The existence of horseshoe dynamics and a coro- 
tation torque is unknown for planets with a horseshoe region width that is a small 
fraction of the turbulent eddy size. In this case, it is possible that turbulence acts 
more as a source of random advection of vortensity through the horseshoe region, 
rather than diffusion. 



4 Migration of gap-opening planets: Type II and III migration 

The disc response to a low-mass planet has been studied in details in §[3] where we 
have focused on the two components of the type I migration torque. The aim of this 
section is to examine the range of planet masses that is relevant to type I migration 
in 



4. 1 and to give a concise description of planet-disc interactions for planets 



that are massive enough to significantly perturb the disc's mass distribution (§ 4.2 
and|431l. 



4.1 Shock formation and gap-opening criterion 

The wakes generated by a planet in a disc carry angular momentum as they prop- 
agate away from the planet. This angular momentum is eventually deposited in the 
disc through some wave damping processes, which leads to redistributing the disc 
mass. An efficient wave damping mechanism relies on the non-linear wave evolu- 
tion of the wakes into shocks |33|. The (negative) angular momentum deposited by 
the inner wake decreases the semi-major axis of the fluid elements in the disc region 
inside the planet's orbit (the inner disc). Similarly, the (positive) angular momentum 
deposited by the outer wake increases the semi-major axis of the fluid elements in 
the outer disc. 

The distance ds from the planet where the planet-generated wakes become shocks 
is given by Il33l 1251 : 

where y is the gas adiabatic index, and a denotes the planet's semi-major axis. 
As the magnitude of the one-sided Lindblad torque peaks at ^ 4H{a)/3 from the 
planet's orbit, a linear description of the differential Lindblad torque thus fails when 
\ds\ ~ 4H{a)/3. This condition can be recast as ^ ^ 0.3h^ for 7 = 5/3. When 
\ds\ < 2i/(a)/3, wakes turn into shocks within their excitation region. Fluid ele- 
ments just outside the planet's horseshoe region are pushed away from the planet 
orbit after crossing the wakes, which directly affects the planet's coorbital region 
by inducing asymmetric U-turns l,66J . Horseshoe fluid elements therefore get pro- 
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gressively repelled from the planet orbit after each U-turn, and the planet slowly 
depletes its coorbital region. The equilibrium structure (width, depth) of the annular 
gap the planet forms around its orbit is determined by a balance between gravity, 
viscous and pressure torques 1 18 1. 

Shock formation and its damping efficiency are very sensitive to the disc's viscos- 
ity, and gap-opening results from a balance between (i) a planet mass large enough 
to induce shocks where the wake excitation takes place, and (ii) a viscosity small 
enough to maximise the amount of angular momentum deposited by the shocks in 
the planet's immediate vicinity: 

1. The first condition reads \ds\ < 2H{a)/3, which corresponds to ^ > 1.5/i^ for 
7=5/3. It means that the planet's Bondi radius rs = GMp/c^, where the pressure 
distribution is most strongly perturbed by the planet, as well the planet's Hill 
radius become comparable to the local pressure scaleheight. This is known as the 
thermal criterion for gap opening ll53l . 

2. The second condition, known as the viscous criterion, can be expressed as ^ > 
40/^, where M = a^^lp/v is the Reynolds numbed] liSTI [131 . 

The above two conditions for gap-opening have been revisited by ITSl . who provide 
a unified criterion that takes the form 

where we have written the disc's kinematic viscosity v = ayh^a^Q OlOOI . and where 



in Eq. ( 34 1 /z and a^ are to be evaluated at the planet's semi-major axis. An illustra- 



tion of the smallest planet mass opening a gap according to criterion ( 34 1 is shown 



in Fig. 13 where it is clear that the gap-opening mass increases with increasing disc 
viscosity and aspect ratio. Assuming h « 0.05, which may be typical of planet form- 
ing regions, the gap-opening mass is in the Saturn-mass range for regions with low 
turbulent activity (dead zones, with typically a^ ^ 10^^), and is in the Jupiter-mass 
range in regions where tty ^ 10^^. Note that when disc self-gravity is included, the 



gap-opening criterion of Eq. ( 34 1 should involve the effective planet mass, that is 
the sum of the planet and circumplanetary disc masses, rather than the planet mass 
alone. 

Prior to a depletion of their coorbital region due to shock formation at the wake's 
excitation region, planets with increasing mass experience a flow transition in their 
immediate vicinity. The flow passes from a low-mass planet configuration described 
3.2.1 and Fig. [3] to a high-mass configuration, where fluid elements may 



m section 



become trapped inside a circumplanetary disc around the planet. This flow transition 
is accompanied by a rapid increase in the half-width Xj of the horseshoe region, 
from a fluid-dominated regime (where x^ °^ [q/h)^'^) to a gravity-dominated regime 



' Although traditionally dubbed Reynolds number essentially for dimensional considerations, this 
ratio has little to do with the dimensionless ratio that must be considered to assess whether a flow 
is laminar or turbulent. If one regards the planet as an obstacle in the sheared Keplerian flow, it 
would be more appropriate to consider as a characteristic scale the size of its Roche lobe or ~ Xs, 
and as a characteristic velocity 2|A|xs. 
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Fig. 13 Minimum planet-to- 
primary mass ratio leading 
to gap-opening as a function 
of the disc's alpha viscous 
parameter (x-axis) and aspect 
ratio (y-axis) at the planet's 
semi-major axis. This mini- 
mum mass is calculated nu- 
merically using criterion l|34|. 
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(where Xs ~ Rh °= q^'^) fTO^. This rapid increase yields a significant increase in the 
corotation torque, as the latter scales as x^. This effect is found to be most significant 
for planet-to-primary mass ratios q ~ 0.6h^ ifTOl . which corresponds to 20 Earth- 
mass planets in h — 0.05 discs. It may contribute to further slowing down, or even 
reversing the migration of growing planets before they carve a gap around their orbit 

MM. 



4.2 Partial gap-opening: type III migration in massive discs 



So far, we have addressed the properties of planet migration through a direct analysis 
of the tidal torque, the latter being directly proportional to the migration rate, see 
Eq. (jTl). This approach is valid for low-mass planets that do not open a gap, for 
which migration has a negligible feedback on the tidal torque (note that a weak, 
negative feedback slightly decreases the magnitude of the entropy-related horseshoe 
drag [68 1). Nevertheless, migrating planets that open a partial gap around their orbit 
experience an additional corotation torque due to fluid elements flowing across the 
horseshoe region |72|. If for instance the planet migrates inwards, fluid elements 
circulating near the inner separatrix of the horseshoe region enter the horseshoe 
region, and execute an outward U-turn when they reach the vicinity of the planet. 
Upon completion of the U-tum, these fluid elements leave the horseshoe region as 
the planet keeps migrating, and end up circulating in the outer disc. Consequently, 
the mass distribution within the horseshoe region may become asymmetric, as the 
horseshoe region adopts approximately a trapezoidal shape in the azimuth-radius 
plane ll72l . As a consequence, in the case of an inward migrating planet, there is 
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more mass behind the planet than ahead of it, owing to the partial depletion of the 



asymmetric horseshoe region. This point is illustrated in the left panel of Fig. 14 
Similarly, if the planet migrates outwards, fluid elements circulating near the outer 
separatrix may embark on single inward U-turns across the horseshoe region. 

Assuming steady migration at a moderate rate (this point will be clarified below), 
the additional corotation torque experienced by the planet due to the orbit-crossing 
flow is, to lowest order in Xsja, 

i^ioss = 1%aaLs X ABaxsi (35) 

where Z^ is the surface density at the inner (outer) horseshoe separatrix for a planet 



migrating inwards (outwards). The term InadEs in the right-hand side of Eq. (35 i 
is the mass flux across the horseshoe region. The second term (ABaxs) is the amount 
of specific angular momentum that a fluid element near a horseshoe separatrix ex- 
changes with the planet when performing a horseshoe U-turn. Note that the above 
expression for i^ross assumes that all circulating fluid elements entering the coor- 
bital region embark on horseshoe U-turns, whereas a fraction of them may actually 
become trapped inside the planet's circumplanetary disc. Since f^ross is proportional 
to, and has same sign as a, migration may become a runaway process. We now 
discuss under which circumstances a runaway may happen. 

The planet and its coorbital material (which encompasses the horseshoe region, 
with mass Mhs, and the circumplanetary disc, with mass M^pd) migrate at the same 
drift rate, a, which we assume to be constant. The rate of angular momentum change 
of the planet and its coorbital region includes (i) the above contribution i^ross to the 
corotation torque, and (ii) the tidal torque that, for planets opening a partial gap, 
essentially reduces to the differential Lindblad torque Ilr- 

2Bad{Mp +Mcpd +Mhs) = 2nadE, x 4Bflx, +rLR. (36) 



Eq. ( 36 1 can be written as 

IBadMp == 2Bad 5m + Flr , (37) 

where Mp — Mp +Mcpd corresponds to an effective planet mass, and where 5m = 
47rflXiZs — Mhs is called the coorbital mass deficit [72 1. It represents the difference 
between (i) the mass the horseshoe region would have if it had a uniform surface 
density equal to that of the separatrix-crossing flow, and (ii) the actual horseshoe 
region mass. 

The migration rate of a partial gap-opening planet, given by Eq. ( [37] i, can be 



described as a feedback loop |66|. This is illustrated in the right panel of Fig. 14 
where the loop input is the differential Lindblad torque, and its output is the mi- 
gration rate. When 5m < Mp, the feedback loop remains stable. The drift rate in 
this case is not strictly a type I nor a type II migration rate. It is rather a type I rate 
enhanced by coorbital effects. No special name has been assigned to this kind of 
migrating regime. When 5m > Mp, the feedback loop gets unstable, and migration 
enters a runaway regime, which can be either inward or outward. The drift rate as 
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Fig. 14 Left: illustration of the flow asymmetry ahead of and behind a Satum-mass planet under- 
going rapid inward runaway migration. Right: type III planetary migration seen as a feedback loop. 
The latter remains stable if the open-loop transfer function £/ x £S < 1, or Sm < Mp. From 1661 . 



a function of the disc mass undergoes a bifurcation ||66l. and this regime is called 
runaway type III migration |7Z 94, 95 1. 

Runaway migration is based on the planet's ability to build up a coorbital mass 
deficit by opening a gap. It does not apply to low-mass planets, for which 5m <C Mp- 
It does not apply to high-mass planets neither, which open a wide, deep gap, so 
that the surface density of the separatrix-crossing flow is too small to produce a 
significant mass deficit. It rather concerns intermediate-mass planets, marginally 
satisfying the gap-opening criterion in Eq. ( [34| , in massive discs (the larger the disc 
mass, the larger the density of the orbit-crossing flow). Its occurrence is illustrated 
in Fig. 15 for a disc with aspect ratio h = 5% and alpha viscosity ay = 4 x 10^^, 
where we see that runaway migration may be particularly relevant to Saturn-mass 
planets in massive discs (with a Toomre-Q parameter at the planet's orbital radius 
typically less than about 10). Bear in mind, however, that the occurrence for runaway 
migration is sensitive to the values of h and ay, since they affect the planet's ability 
to open a partial gap. Also, note that the criterion for runaway migration features the 
effective planet mass Mp, sum of the planet mass and the circumplanetary disc mass. 
The occurrence for runaway migration is therefore sensitive to the mass distribution 
inside the circumplanetary disc, which may be significantly affected by the assumed 
physical modelling, e.g., whether gas accretion on the planet is taken into account 
124 1, the inclusion of self-gravity [112 J, or the treatment for the gas thermodynamics 
|f93J . It may also be affected by grid resolution effects in numerical simulations 



The simple model described above helps understand the condition for migra- 
tion to enter a runaway regime. However, since it assumes steady migration (that 
is, constant a), this model is no longer valid when migration actually enters the 
runaway regime, where the migration rate increases exponentially over a time com- 
parable to the horseshoe libration period. A more general approach can be found 
in ll72l l66l l94l . As long as the orbital separation by which the planet migrates over 
a Ubration period remains smaller than the radial width of the planet's horseshoe 
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Fig. 15 Occurrence for type I, II and III (runaway) migrations with varying the planet-to-primary 
mass ratio (bottom x— axis) and the disc-to-primary mass ratio at the planet location (left y— axis). 
The disc's aspect ratio is h = 0.05 and its alpha viscosity is tty = 4 x 10^^^. The right >■— axis 
shows the Toomre Q-parameter at the planet location. The upper part of the plot is limited by the 
gravitational instability limit (dashed line). From 1 72] , 



region, i^ross remains approximately proportional to a (slow runaway regime). At 
larger migration rates (fast runaway regime), f^ioss reaches a maximum and slowly 
decreases with increasing a [72| (see also Fig. 16 in |66|). The precise dependence 
of i^ross with a in this fast runaway regime is intrinsically related to the evolution 
of the mass coorbital deficit, and therefore to the planet's migration history. The or- 
bital evolution of planets subject to runaway type III migration is therefore difficult 
to predict. Numerical simulations find that, depending on the resolution of the gas 
flow surrounding the planet, the timescale for inward runaway type III migration 
can be as short as a few 10^ orbits |72, 15 1. 

The sign of f^ross is dictated by the initial drift of the planet. Runaway migration 
can therefore be directed inwards or outwards, depending on the sign of a before the 
runaway takes place. In particular, migration may be directed outwards if, despite 
the coorbital region being partly depleted, a (positive) horseshoe drag remains strong 
enough to counteract the (negative) differential Lindblad torque. Outward runaway 
migration could thus be an attractive mechanism to account for the recent discovery 
of massive planets at large orbital separations (which we will further discuss in 

however show fliat, despite the expected increase in 



5.1 1. Simulations 



the mass of the orbit-crossing flow as the planet moves outwards (for background 
surface density profiles shallower than r^^), the mass coorbital deficit cannot be 
retained indefinitely. The increase in the mass of the circumplanetary disc, and the 
strong distortion of the flow within it at large migration rates, lead the planet to 
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eventually lose its coorbital mass deficit, and the sense of migration is found to 
reverse. 

Type III migration has been recently revisited in low-viscosity discs (a-j < 
a few X 10^^). Depending on the disc mass, the edges of the planet-induced gap may 
be subject to two kinds of instabilities. In low-mass discs, gap edges are unstable to 
vortex-forming modes 1581 l50l |56 1 . They lead to the formation of several vortices 
sliding along the gap edges, which merge and form large-scale vortices. When they 
pass by the planet, these vortices may embark on horseshoe U-turns and exert a large 
corotation torque on the planet, with the consequence that the planet can be scattered 
inwards or outwards 1541 . When the fluid's self-gravity is taken into account, only 
a fraction of the large-scale vortices actually embark on horseshoe U-turns, the rest 
of the vortices keeps on sliding along the gap edges. This provides a periodic, inter- 
mittent corotation torque on the planet. Depending on the relative strengths of the 
vortices embarking on inward and outward U-turns, this mechanism acts much like 
an intermittent type III migration regime. In massive self-gravitating discs (stable 
against the gravitational instability), vortex-forming modes are replaced by global 
edge modes, which excite spiral density waves ll55l . A decreasing radial profile of 
the Toomre-Q parameter favours edge modes at the gap's outer edge. The periodic 
protrusion of edge mode-induced density waves near the gap's outer edge provides 
a periodic source of (positive) corotation torque on the planet, and induces an in- 
termittent type III migration regime. Numerical simulations by |57| show that edge 
modes can sustain outward migration, until the planet leaves its gap. 

We briefly sum up the main results of this paragraph: 

• Migrating planets experience an additional corotation torque due to fluid ele- 
ments flowing across the horseshoe region, and embarking on horseshoe U-turns. 
It is proportional to the planet's migration rate at small migration rates, which 
gives a positive feed back on migration. When the feedback loop diverges, the 
migration type is known as type III migration. 

• The planet and its circumplanetary disc feel an effective corotation torque that 



is proportional to the coorbital mass deficit, defined through Eq. (37i. The oc- 
currence of a runaway feedback (i.e. of type III migration) corresponds to the 
coorbital mass deficit exceeding the mass of the planet and its circumplanetary 
disc. This applies to planets opening a partial gap around their orbit in massive 
protoplanetary discs. 

The orbital evolution of planets undergoing type III migration is sensitive to 
the time evolution of the coorbital mass deficit, which makes it difficult to pre- 
dict. Numerical simulations show that runaway migration operates on very short 
timescales, typically in 100 to 1000 planet orbits. 



4.3 Deep gap-opening: Type II migration 

Planets massive enough to clear their coorbital region and open a deep gap around 



their orbit (see illustration in Fig. 16 1 enter the migration regime called Type II 
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Fig. 16 Gap opened by a 
Jupiter-mass planet orbiting a 
Sun-like star. 




migration. Such planets satisfy the gap-opening criterion given by Eq. ( 34 1. Assum 



ing, for instance, a protoplanetary disc with aspect ratio ~ 5% and alpha viscosity 
ttv ^ 10^^, type II migration typically applies to planets more massive than Jupiter 
orbiting Sun-like stars. Compared to the type I and type III migration regimes de- 
scribed previously, the amplitude of the corotation torque is much reduced due to 
the clearing of the planet's coorbital region, and the differential Lindblad torque 
balances the viscous torque exerted by the disc. The net torque on the planet can 
be written as a fraction Cu of the viscous torque due to the outer disc ifTTl . This 
factor Co features the time-dependent fraction of gas /gas remaining in the planet's 
coorbital region. 

The particular case with /gas going to zero corresponds to what is usually referred 
to as standard type 11 migration regime. Its timescale can be approximated as 

where v denotes the disc's kinematic viscosity, E the surface density of the disc 
perturbed by the planet, and r^ is the location in the outer disc where most of the 
planet's angular momentum is deposited. It can be approximated as the location of 
the outer separatrix of the planet's horseshoe region, which, for gap-opening planets. 



is ^o ~ a + 2.5RH [?0J. The first term on the right-hand side of Eq. (38 1 corresponds 
to the viscous drift timescale at radius Tq, and the second term features the ratio 
of the planet mass to the local disc mass at radius r^. Two migration regimes can 
therefore be distinguished: 

I . Disc-dominated type II migration. When the planet mass is much smaller than the 
local disc mass (by which we refer to the quantity AnE {ro)r^), the planet behaves 
much like a fluid element that the disc causes to drift viscously. The planet's mi- 
gration timescale then matches the disc's viscous drift timescale « 2rl/3v{ro) 
ll52l . In this migration regime, called disc-dominated type II migration, the planet 
remains confined within its gap. Should the planet migrate slightly faster than 
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the disc near its orbit, the increased inner Lindblad torque due to the planet get- 
ting closer to the gap's inner edge would push the planet outward. Conversely, 
should the planet migrate at a slower pace than the (local) disc, the increased 
outer Lindblad torque would push the planet back inward. The timescale for the 
disc-dominated type II migration regime, Tm, can be recast as 

.„,«4.7xl0Vsx(^)"f^y7'^y"V^)''', (39) 
' ^ VlO-2/ \0.Q5j \MqJ \5A\jJ ' ^ ^ 

where a^ and h are to be evaluated at Tq. 

2. Planet-dominated type II migration. When the planet mass becomes comparable 
to, or exceeds the local disc mass, the orbital evolution of a gap-opening planet 
is no longer dictated by the disc alone. The inertia of the planet slows down 
its orbital migration BlOll [39l . and in the limit when the planet mass is large 
compared to the local disc mass, the planet enters the so-called planet-dominated 
type II migration regime, whose timescale Tu p reads 



with Tiid given by Eq. ( 39 1, and where the planet to local disc mass ratio reads 



^P ■.nnf^p\f ^('•o) 



-2 



47rZ(ro)r2~^°^(,M0J 1^150 gcm-2J IsAuJ " ^^^^ 

In self-gravitating discs, the distinction between the planet- and disc-dominated 
type II migration regimes should involve the comparison between the local disc 
mass and the effective planet mass Mp, that is the sum of the planet and circumplan- 
etary disc masses. Related to this point, we comment that the planet and its circum- 
planetary disc migrate at the same drift rate. When its self-gravity is included, the 
protoplanetary disc torques both the planet and the circumplanetary disc. However, 
if self-gravity is discarded, as is usually the case in numerical simulations, the pro- 
toplanetary disc can only torque the planet, and the circumplanetary disc remains a 
passive spectator of the migration. In this case, the planet must exert an additional 
effort to maintain the planet and circumplanetary disc joint migration. Put another 
way, the circumplanetary disc artificially slows down migration when self-gravity 
is discarded. To avoid this artificial slowdown, ifTSl showed that, in simulations dis- 
carding self-gravity, the calculation of the torque on the planet must exclude the ma- 
terial inside the circumplanetary disc. In addition, migration rates with and without 
self-gravity can be in close agreement, provided that the mass of the circumplane- 
tary disc is added to that of the planet when calculating the gravitational potential 
felt by the protoplanetary disc. 

In the early stages of their formation and orbital evolution, most massive gap- 
opening planets should be subject to disc-dominated type II migration, and migrate 
on a timescale comparable to the disc's viscous timescale. Note from Eq. ([38]l that 
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this corresponds to the shortest migration timescale a gap-opening planet can get. 
Depletion of the protoplanetary disc, or substantial migration towards the central 
object, should, however, slow down migration as the planet's inertial mass becomes 



comparable to the local disc mass. It is nonetheless interesting to note from Eq. ( 39 1 
that in the early stages of the disc evolution, type II migration can be relatively fast 
in the disc's turbulent parts. This may make difficult the maintenance of massive 
planets at reasonably large orbital separations from their host star. Additional mech- 
anisms, like the effect of stellar irradiation on the disc's density and temperature 
profiles near the planet's orbit (formation of shadow regions near the gap's inner 
edge, and irradiated 'puffed up' regions near the gap's outer edge) |40| could help 
slow down type II migration. 

Gap formation and type II migration are intimately related to the disc viscosity 
in laminar viscous disc models, and a few studies have investigated their proper- 
ties in MHD turbulent discs lISTl IIIOI lOTI . These studies have considered three- 
dimensional magnetised disc models, where vertical stratification and non-ideal 
MHD effects were discarded for simplicity. They found that the structure of the 
annular gap opened by a massive planet in fully MHD turbulent discs is essentially 
in line with the predictions of viscous disc models with a similar alpha viscous pa- 
rameter near the planet location. Gaps in turbulent discs tend, however, to be wider 
than in viscous discs llSTlfTTO I. Some other differences arise between turbulent and 
viscous disc models, particularly in the vicinity of the planet, where magnetic field 
lines are compressed and ordered at the location of the wakes and the circumplan- 
etary disc. The connection between the circumplanetary and protoplanetary discs 
through magnetic field lines can cause magnetic braking of the circumplanetary 
material |91 1, which may help increase gas accretion onto the planet 1 110, 91 1. 

Before leaving this section, we comment that the overall properties of planet-disc 
interactions with gap-opening planets remain essentially unchanged when taking the 
disc's vertical stratification into account, and therefore the two-dimensional approx- 
imation is valid. Nonetheless, different structures in the flow circulating around the 
planet in two- and three-dimensions, and the related accretion rate onto the planet, 
may affect the planet's migration rate (see e.g., |2()|). 

We summarise the results described in this section: 

• Planets massive enough to open a wide and deep annular gap around their orbit 
are subject to Type II migration. 

• When the local disc mass (roughly speaking, the mass interior to the planet or- 
bit) remains large compared to the planet mass, the planet migration timescale 
corresponds to the viscous drift timescale (disc-dominated type II migration). 

• When the planet mass becomes comparable to, or exceeds the local disc mass, 
migration is slowed down by the planet's inertia (transition to planet-dominated 
type II migration). 
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5 Planet migration theories and exoplanets' observed diversity 

The properties of planet-disc interactions have been examined in details in the pre- 
vious sections, with a particular emphasis on the expected migration rate for planets 
with different masses. The migration rate is intimately related to the disc's physical 
properties (e.g., mass, sound speed, cooling properties, turbulent stresses) near the 
planet location, which underlines that the modelling of protoplanetary discs plays 
as much of an important role as planet-disc interaction theories in predicting the 
evolution of planetary systems. We continue our exposition with a brief discussion 
of several aspects of planet-disc interactions which could account for the observed 
diversity of (extrasolar) planets. 



5.1 Massive planets at large orbital separations 

Amongst the recent discoveries of exoplanets, particularly exciting is the observa- 
tion by direct imaging of about 10 massive exoplanets located at separations ranging 
from 10 to 200 AU from their host star (e.g., [.61..47J ). Most of these planets are so 
far observed to be the only planetary companions of their host star. Yet, it is possible 
that their present location results from a scattering event with another massive com- 
panion on a shorter-period orbit. A remarkable exception is the HR 8799 planetary 
system. It comprises four planets with masses evaluated in the range [7 — 10] Jupiter 
masses, and estimated separations of 14, 24, 38 and 68 AU fST]. The planets are 
close to being in mutual mean-motion resonances, and it seems likely that planet- 
disc interactions could have played a major role in shaping this planetary system. 
We discuss below the relevance of planet-disc interactions to account for massive 
planets at large orbital separations. 



5.1.1 Outward migration of a pair of massive resonant planets 

In the standard core -accretion scenario for planet formation, it is difficult to form 
Jupiter-like planets in isolation further than ^10 AU from a Sun-like star 1.971 137]| . 
As we have seen in § B] planets in the Jupiter-mass range orbiting Solar-type stars 
are expected to open an annular gap around their orbit. If a partial gap is opened, 
outward runaway type III migration could occur under some circumstances, but 
as we have discussed in § |4.2| numerical simulations indicate that it is difficult to 
sustain this outward migration in the long term. If the planet opens a deep gap, 
inward type II migration is expected. It is therefore unlikely that a single massive 
planet formed through the core-accretion scenario within ^ 10 AU of its host star 
could migrate to several tens of AUs. 

A notable exception to this generally expected trend has been recently proposed 
by llT6l . based on a migration mechanism originally studied by 163)1 . This mecha- 
nism relies on the joint migration of a pair of resonant massive planets embedded 
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in a common gap. In this mechanism, the innermost planet is massive enough to 
open a deep gap and migrate inwards on a timescale comparable to that of type 
II migration. The outermost, less massive planet migrates inwards at a larger pace 
while carving a partial gap around its orbit. If both planets open overlapping gaps, 
and maintain a mean-motion resonance between their orbits, their joint migration 
could proceed outwards. The global picture is the following: as the inner planet is 
more massive, the torque it experiences from the inner disc (inner Lindblad torque) 
is larger than the (absolute value of the) torque the outer planet experiences from the 
outer disc (outer Lindblad torque). To maintain joint outward migration in the long 
term, the fluid elements outside the common gap must be funnelled to the inner disc 
by embarking on horseshoe trajectories. Otherwise, material would pile up at the 
outer edge of the common gap, much like a snow-plough, and the torque balance as 
well as the sense of migration would eventually reverse. An illustration of the joint 



outward migration mechanism is shown in Fig. 17 




Fig. 17 Illustration of the joint outward migration of a pair of resonant massive planets. The left 
panel shows the evolution of the disc's surface density perturbed by a Jupiter-mass planet (inner 
planet) and a Saturn-mass planet (outer one). After an episode of rapid convergent migration (top- 
left quadrant) resulting in their capture into mean-motion resonance, planets open overlapping 
gaps (top-right quadrant), which leads to their joint outward migration (lower quadrants). The right 
panels illustrate the outcome of the same mechanism applied to an inner 3-Jupiter mass planet, and 
an outer 2- Jupiter-mass planet orbiting a 2M3 mass star (taken from [16;|). The time evolution of 
the planets semi-major axis is shown in the top-right panel, and that of their 2 : 1 critical resonance 
angle is in the bottom-right panel. 



The migration reversal described above requires an asymmetric density profile 
within the common gap. It is thus sensitive to the disc's aspect ratio and viscosity, 
which enter the gap-opening criterion. It is also sensitive to the mass ratio of the 
two planets. If the outer-to-inner planets mass ratio is too small, the density contrast 
within the common gap will be too large to affect the evolution of the innermost 
planet (the gas density near the outer planet's orbit remains too large to signifi- 
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candy decrease the outer Lindblad torque acting on the inner planet). Conversely, 
if the outer-to-inner planets mass ratio is too large, the Lindblad torques imbalance 
will favour joint inward migration. By changing the planets mass ratio during joint 
migration, gas accretion onto the planets could affect the possibility of sustaining 
outward migration in the long term. This issue requires further investigation, and 
accurate modelling of the gas accretion processes onto Saturn sized planets. 

We also mention that the joint outward migration scenario has been recently dis- 



cussed in the context of the Solar System II105I . Inward migration of Jupiter in 
the primordial Solar nebula down to « 1.5 AU, followed by joint outward migration 
with Saturn to the current location of both planets (the "Grand tack") would truncate 
the disc of planetesimals interior to Jupiter's orbit at about 1 AU. The subsequent 
formation of the terrestrial planets is found to occur with the correct mass ratio be- 
tween Earth and Mars, and would also account for the compositional structure of 
the asteroid belt Ml 051. 



5.1.2 Migration of planets formed by gravitational instability 

An alternative to the core-accretion formation scenario involves the fragmentation 
of massive protoplanetary discs into clumps through the gravitational instability. 
Gravitational instability (GI) may typically occur at separations larger than 30 to 
50 AU from a central (Sun-like) star, if the Toomre-Q parameter approaches unity 
and the disc's cooling timescale becomes of order the dynamical timescale (e.g., 
||98ll29l ). While several massive planets could form by fragmentation of a massive 
disc at several tens of AUs from their star, they are unlikely to stay in place. The 
tidal interaction with the gravito-turbulent disc they are embedded in should rapidly 
bring planets formed by GI to the disc's inner regions llll l75l[TT3l . in a timescale 
comparable to that of type I migration [ 11 J . The orbital evolution of a single Jupiter- 
mass planet embedded in a gravito-turbulent disc (where the planet is supposed to 



have formed by GI) is illustrated in Fig. 18 where we see that the planet migrates 



from 100 to 20 AU in typically less than 10 yrs. 

Investigation is under way to determine the evolution of planets formed by GI 
when they reach the inner parts of protoplanetary discs. The latter should be too hot 
to be gravitationally unstable, and other sources of turbulence, such as the magne- 
torotational instability, could prevail, changing the background disc profiles as well 
as the amount of turbulence. It is thus possible that the rapid type I migration of 
planets formed by GI slows down in the disc inner parts and results in the formation 
of a gap. Gap-opening may also occur if significant gas accretion occurs during the 
planets fast inward migration 11131 . Planet-planet interactions, which may result in 
scattering events, mergers or captures in resonance, should also play a prominent 
role in shaping planetary systems formed by GI. 
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Fig. 18 Jupiter-mass planet embedded in a gravito-turbulent disc. After setting up a quasi 
steady-state gravito-turbulent disc with (gravito-turbulent) shock heating balancing disc cooling 
(parametrised here by a simple /3-cooling function, see 1 11 1), simulations were restarted with in- 
serting a Jupiter-mass planet at 100 AU. The left panel shows the disc's surface density three orbits 
after restart, the planet being located at x ~ — 50 AU, j ~ 60 AU. The right panel displays the 
time evolution of the planet's orbital separation in 8 different restart simulations with varying the 
azimuth of the planet prior to its insertion in the disc. Taken from 1111 . 



5.2 Planet population syntheses 



Planetary astrophysics is undergoing an epoch of explosive growth, driven by the ob- 
servational discoveries of more than 750 exoplanets over the past two decades. Out- 
standing progress in detection techniques have uncovered planetary systems very 
different from ours. Since the discovery of the first Hot Jupiter |73 1, radial velocity 
surveys have made possible the detection of Earth-like planets, some in the habitable 
zone of their star ||92l . Transit space missions CoRoT and KEPLER are digging out 
hundreds of close-in extrasolar planets, some in exotic environments (like Kepler- 
16 b, the first circumbinary exoplanet discovered |26|). Direct imaging has revealed 
the existence of massive giant planets located at several tens of AU from their star 
Such diversity provides an exciting opportunity to test our theories for the forma- 
tion and evolution of planetary systems. By coupling theoretical models of planet 
formation and migration, and of disc evolution, planet population syntheses estimate 
the statistical distribution of exoplanets according to their mass, semi-major axis, 
and eccentricity, which they compare to observed distributions [37, 38, 76 99, 35]. 
At the moment, models of planet population syntheses are not able to reproduce 
the statistical properties of extrasolar planets. For instance, they predict a deficit of 
super-Earths and Neptune-like planets with orbital periods less than 50 days, while 
observations have revealed a significant number of exoplanets in this range of mass 
and period |361. The origin for this discrepancy can be found in uncertain prescrip- 
tions for the minimum core mass for the onset of gas accretion ||35]| . as well as in 
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the modeling of type I migration. The difficulty raised by the excessively rapid in- 
ward type I migration, predicted by the long-time reference torque formula by 11021 . 
was circumvented by introducing a reduction factor in front of this torque formula. 
Population syntheses models tried to constrain this factor to reproduce the statisti- 
cal properties of detected exoplanets. This reduction factor was found to range from 
0.01 to 0. 1 . We note however that the introduction of this reduction factor to provide 
planetary population synthesis in agreement with the statistics of detected extraso- 
lar planets is ad hoc, and that there is no reason to expect that the type I migration 
drift rates are systematically overestimated in theoretical studies by a factor 10 to 
100. Rather, as we emphasized throughout this manuscript, type I migration is very 
sensitive to the disc's density and temperature profiles near the planet orbit. Large 
slopes of mass density and/or temperature, over a limited radial range, may reverse 
the tidal torque exerted on the planet. This, in turn, may create "planetary traps" at 
the points where the tidal torque cancels out (much like what was contemplated by 
ITTTI for the case of a positive surface density gradient), which can stop incoming 
protoplanets, depending on their mass. The number and location of these traps may 
vary as the disc evolves. This view of type I migrating objects subject to several 
traps on their way to the star ||60ll34ll sounds more compatible with the state of mi- 
gration theories than an ad hoc reduction factor. This has motivated recent works 
to produce accurate, yet simple formulae for type I migration ll69l [85l l67l . These 
formulae include a description of the corotation torque in discs with arbitrary vis- 
cosity and thermal diffusion, and corrections to the Lindblad torque for discs with 
non power-law profiles. Their incorporation into models of planet population syn- 
thesis will hopefully provide a better comprehension of the diversity of observed 
exoplanets. 



6 Conclusions 

We have reviewed the recent progress made in understanding planet-disc interac- 
tions, and the resulting planets' orbital migration. We have particularly focused on 
the migration of growing protoplanets (type I migration), which has been the subject 
of intensive investigation over the past five years. Being for a while the second-place 
actor of planet migration theories, the corotation torque has been shown to play a 
prominent role in realistic protoplanetary discs, where it can slow down, stall, or 
reverse type I migration. This review is especially aimed at giving a comprehensive, 
detailed description of the mechanisms responsible for the corotation torque. The 
type II and type III migration regimes for gap-opening planets are also reviewed and 
discussed in the context of observed exoplanets. Being aimed at migration of plan- 
ets on circular orbits, this review has set aside interesting recent developments on 
the tidal interactions of eccentric or inclined planets with their discs. We have also 
focused essentially on the mechanisms that drive the migration of a single planet 
in a disc, and we have therefore excluded most results about the migration of sev- 
eral planets. For a recent review covering these topics, the reader is referred to Ii45]| . 
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This list of restrictions of the present review stresses that the research on planet- 
disc interactions is a very active branch of planet formation, with a growing body 
of avenues. We finally reiterate the plea made in the introduction: planetary migra- 
tion is not overrated. The tremendous value of each of the tidal torque components 
exerted on a given planet, associated to the great sensitivity of these torques to the 
underlying disc structure, appeals for a detailed knowledge of the properties of pro- 
toplanetary discs, and significant efforts toward an accurate determination of each 
torque component. This also reasserts tidal interactions as a prominent process in 
shaping forming planetary systems. 
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